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IV 


ABSTRACT 


The  Naval  Postgraduate  School’s  Ultraviolet  Imaging  Spectrometer  (NUVIS)  is 
a  hyperspectral  sensor  with  a  spectral  response  of  300  to  375  nanometers.  This 
thesis  research  concentrates  on  the  operation  and  calibration  of  NUVIS  in  the 
detection  of  effluent  sulfur  dioxide  (SO2)  plumes.  NUVIS  is  capable  of  detecting 
and  quantifying  SO2  emissions  in  the  form  of  effluent  smokestack  plumes  by 
exploiting  S02’s  unique  UV  absorption  signature.  Laboratory  comparison  UV 
spectra  of  SO2  were  recorded  and  used  to  calculate  curves  of  growth  for  four 
different  SO2  spectral  features.  Laboratory  results  were  employed  to  analyze  field 
data  taken  of  a  coal-buming  power  plant.  Analysis  of  this  plume  data  yielded  a 
mean  plume  SO2  mixing  ratio  of  365  ±  200  ppm,  in  agreement  with  the  in  situ  stack 
value  of  400  ppm.  Further  assessment  of  NUVIS  indicates  that  its  lower  limit  for 
SO2  detection  in  typical  field  applications  is  approximately  70  ppm. 
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INTRODUCTION 


A.  MOTIVATION  AIVD  OBJECTIVES 

Remote  sensing  instruments  are  critical  to  the  areas  of  technical  intelligence,  the 
remote  sensing  of  hazards,  and  support  to  military  operations.  In  an  effort  to  enhance  the 
capabilities  available  to  the  remote  sensing  community  as  well  as  to  the  military,  we  have 
pursued  a  little  used  portion  of  the  electromagnetic  spectrum.  Over  the  last  three  years, 
we  have  designed  and  constructed  an  ultraviolet  hyperspectral  imager.  Our  motivation 
came,  in  part,  from  the  success  of  the  school’s  Middle  Ultraviolet  SpecTrograph  for 
Analysis  of  Nitrogen  Gases  (MUSTANG),  a  rocket-home  atmospheric  research 
instrument.  This  resulted  in.  the  development  of  the  NFS  UltraViolet  Imaging 
Spectrometer  (NUVIS),  a  hyperspectral  imager  that  operates  at  near-ultraviolet 
wavelengths.  NUVIS  is  one  of  the  first  such  instruments  of  its  kind:  a  portable,  ground- 
based,  hyperspectral  imager  that  operates  at  ultraviolet  (UV)  wavelengths:  nominally  300 
to  375  nanometers  (nm). 

From  the  onset,  NUVIS  designers  recognized  that  ultraviolet  remote  sensing 
could  be  used  as  an  important  tool  in  effluent  analysis.  Through  the  efforts  of  Dr.  David 
Cleary,  LT  Todd  Hooks  and  Major  Andrew  MacMannis,  NUVIS  was  designed  at  the 
school,  and  fabrication  of  the  instrument  was  completed  in  the  fall  of  1997.  The  next  step 
was  to  focus  research  on  the  instrument’s  ability  to  analyze  and  characterize  a  chemical 
effluent  plume.  Such  capabilities  are  of  interest  to  both  the  civilian  and  military 
communities. 
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As  initial  testing  was  about  to  get  underway,  a  major  question  arose:  Could 
gaseous  effluents  emitted  from  smoke  stacks  be  detected,  imaged,  identified,  and 
quantified  using  hyperspectral  remote  sensing  techniques  at  UV  wavelengths?  This 
thesis  research  set  out  to  answer  that  question  and  to  further  define  the  instrument’s 
capability  to  perform  such  a  mission.  In  support  of  these  efforts,  the  trace  gas,  sulfur 
dioxide  (SO2),  was  targeted  for  study. 

Numerous  studies  have  focused  on  the  topic  of  remote  sensing  of  trace  gases 
using  spectroscopic  monitoring  techniques.  Supporting  literature  indicates  that 
considerable  attention  has  been  placed  on  the  study  and  techniques  of  infrared  remote 
sensing.  Until  now,  remote  sensing  at  ultraviolet  wavelengths  has  received  compartively 
little  attention.  This  study  sets  out  to  challenge  such  thinking  by  exploring  new  ways  to 
exploit  the  ultraviolet  region  of  the  electromagnetic  spectrum  utilizing  the  recently 
developed  technique  of  hyperspectral  imaging. 

Although  ultraviolet  wavelengths  suffer  more  from  atmospheric  scattering  and 
ozone  absorption  than  do  others,  this  study  will  show  that  NUVIS  has  already 
demonstrated  dual-use  imaging  and  spectroscopic  capabilities  and  is  capable  of  detecting, 
at  short  range,  the  effluent  trace  gas  SO2.  With  this  initial  success,  it  became  necessary  to 
determine  whether  effluent  SO2  could  reliably  be  identified,  measured  and  quantized  over 
longer  distances  using  UV  remote  sensing.  An  approach  was  taken  such  that  laboratory 
experiments,  calibrations,  and  operational  field  testing  would  lead  to  an  evaluation  of 
NUVIS ’s  capabilities  to  detect  the  trace  gas  and  the  utility  of  such  technology. 
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This  thesis  project  focuses  on  calibrating  the  spectrometer  for  the  detection  of 
gaseous  sulfur  dioxide.  Laboratory  experiments  were  conducted  at  known  concentrations 
of  SO2  to  measure  changes  in  its  absorption  spectrum  as  a  function  of  gas  concentration. 
SO2  is  a  triatomic  molecule  that  has  a  set  of  strong  absorption  signatures  between  300 
and  325  nm.  The  detailed  spectrum  across  this  spectral  region  exhibits  a  number  of 
identifiable  absorption  lines.  For  each  observed  line,  an  associated  measure  of  its 
spectral  absorption,  called  its  equivalent  width,  may  be  easily  measured.  Calibration 
curves  of  the  lines’  equivalent  widths  versus  gas  concentration,  called  curves  of  growth, 
can  then  be  generated. 

During  this  research,  six  sets  of  calibration  data  were  analyzed.  From  this 
analysis,  four  individual  absorption  lines  were  identified  for  equivalent  width 
measurements.  Four  curves  of  growth  were  successfully  generated  from  these  laboratory 
experiments.  These  calibration  curves  were  then  applied  to  field  data  in  an  initial 
attempt  to  quantify  effluent  SO2  concentrations. 

It  should  be  noted  that  my  involvement  with  NUVIS  came  at  a  time  when  nearly 
all  of  the  instrument’s  initial  designers  had  either  moved  on  from  the  school  or  were  in 
the  process  of  leaving.  Because  of  this,  a  period  existed  where  Dr.  Cleary  had  done  a 
considerable  amount  of  work  with  the  instrument  himself,  and  the  majority  of  this  work 
had  gone  unpublished.  Because  Hooks  and  MacMannis  were  unable  to  complete 
laboratory  and  field  testing  of  the  instrument,  a  detailed  operating  procedure  had  not  yet 
been  developed.  This  thesis  will  not  only  present  the  work  involved  with  generating  the 


SO2  curves  of  growth  but  will  also  outline  the  basic  operating  procedures  for  the 
instrument.  Such  documentation  had  not  existed  prior  to  this  study. 

In  an  attempt  to  maintain  continuity  with  the  work  done  by  Hooks  and 
MacMannis,  this  thesis  also  includes  certain  portions  of  Dr.  Cleary’s  unpublished  work 
as  backgrormd  material.  This  is  in  part  due  to  the  fact  that  the  current  NUVIS  system  has 
undergone  certain  improvements  from  its  original  configuration.  Where  appropriate,  the 
inclusion  of  such  material  is  clearly  indicated. 

In  1998,  ^ound-based  field  experiments  were  conducted  and  resulted  in 
preliminary  spectroscopy  and  spectral  imagery  of  effluent  S02-producing  targets.  This 
thesis  research,  the  accompanying  calibration,  a  follow-on  field  experiment,  and  the 
resultant  data  analysis  have  all  been  accomplished  in  support  of  chemical  plume 
detection,  analysis  and  characterization  with  respect  to  sulfur  dioxide.  The  work 
involved  in  this  study  is  a  follow-on  to  the  initial  experimental  research  initiated  in  1998. 

B.  OUTLINE 

This  thesis  is  divided  into  six  chapters  and  three  appendices.  Chapter  11  provides 
the  necessary  background  material  on  the  topics  of  remote  sensing  and  hyperspectral 
imaging,  ultraviolet  remote  sensing,  sulfur  dioxide  and  absorption  spectroscopy.  The 
NIXVIS  instrument  will  be  presented  in  Chapter  III.  Chapter  IV  details  the  laboratory 
calibration  of  NUVIS  and  includes  resultant  curves  of  growth,  one  for  each  of  the  four 
SO2  absorption  lines  selected  for  analysis.  Field  data  is  briefly  analyzed  by  a  simple 
application  of  these  laboratory  calibrations  developed  in  Chapter  IV.  This  process  and 
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the  accompanying  results  are  presented  in  Chapter  V.  Conclusions  and  recommendations 
are  contained  in  Chapter  VI.  Appendix  A  provides  a  recommended  checklist  for  NUVIS 
field  operation.  Appendix  B  briefly  details  the  system  operating  procedure.  Appendix  C 
provides  additional  data  that  should  be  considered  in  future  reasearch. 
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II.  BACKGROUND 


A.  REMOTE  SENSING  AND  SPECTRAL  IMAGING 

Remote  sensing  is  a  generic  term  used  to  describe  all  of  the  measurements  made 
at  a  distance  from  an  object  to  be  observed.  It  involves  collecting  information  about  an 
object  without  being  in  direct  contact  with  that  object.  Typically  the  information  being 
collected  is  carried  by  electromagnetic  waves,  although  other  means  are  sometimes  used. 
If  the  source  of  these  waves  is  a  man-made  emitter,  then  the  process  is  dubbed  active 
remote  sensing.  If  the  source  of  these  waves  is  natural,  the  process  is  known  as  passive 
remote  sensing  (Lenoble,  1993).  It  is  passive  remote  sensing  at  near-ultraviolet 
wavelengths  (200-400nm)  with  which  NUVIS  is  most  concerned.  As  a  quick  reference, 
Table  2.1  lists  some  typical  wavelength,  frequency  and  wavenumber  ranges  for  various 
regions  of  the  electromagnetic  spectrum. 


Name 

Wavelengths 

Frequencies 

Wavenumbers 

(nm) 

(Hz) 

(cm-') 

Radio 

io-®-o.i 

Microwave 

-3  X  10''-3  X  10'" 

0.1-1 

Infrared 

-3  X  10'"-3  X  10''’ 

1-10' 

Visible 

700-400 

-(4.3-7.1)x  10''' 

(1. 4-2.4)  X  10' 

Ultraviolet 

-400-50 

-7.5  X  10'^-6.0  X  10" 

■ nil— 11 

XRay 

-50-0.1 

-(0.6-300)  X  10’® 

■nEBrimsriaii^^H 

Gamma  Ray 

<0.1 

To  3  X  10" 

>  10* 

Table  2.  1 .  Typical  wavelengths,  frequencies  and  wavenumbers  for  various  regions  of  the 

electromagnetic  spectrum. 
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Over  the  last  decade,  a  number  of  instrument  platforms  have  been  developed  for 
the  express  purpose  of  optical  remote  sensing.  These  platforms  have  been  employed  in  a 
number  of  field  operations,  ranging  from  space-based  and  air-based  missions  to  ground- 
based  applications.  Spectral  imagers  are  one  such  instrument  family,  and  they  have 
quickly  become  v^orkhorses  of  the  remote  sensing  community.  Spectral  imagers  measure 
both  the  2-D  image  produced  by  a  scene  and  the  spectral  content  of  each  picture  element 
(pixel)  of  the  image.  These  spectral  imagers  are  classified  by  the  width  of  the  spectral 
bands  they  record.  Figure  2.1  provides  a  brief  overview  of  such  classifications. 

The  motivation  behind  the  development  of  NUVIS  was  a  perceived  need  to 
explore  the  near-UV  wavelength  region  with  an  operational  hyperspectral  imager.  Such 
an  instrument  would  be  capable  of  collecting  data  from  tens  to  hundreds  of  bands  at  a 
much  smaller  bandwidth  resolution  than  that  of  a  conventional  multispectral  sensor. 


WAVELENGTH  (ran) 


MULTISPECTRAL 
Wide  Bandwidth 
IVIoderate  Spectral  Resolution 


HVPERSPECTRAL 
Narrow  Bandwidth 
High  Spectral  Resolution 


ULTRASPECRAL 
'  Very  Narrow  Bandwidth 
Very  High  Spectral  Resolution 


400  700 


1000  2000 


4000 


6000  14000 


ftmiiiijiiiiii  . 


l^cctral  identification  of  major  features, 
i;e.,  trees,  grass,  roads 


liil 


Spectral  discrimination  by  species , 
materials,  and  environmental  conditions 


Identification  and  discrimination  of  subtle 
ilspectral  details  of  materials,  v^ors  and  aerosols 


Figure  2.  1.  Spectral  imager  classifications.  From  Multispectral  Users  Guide,  1995. 
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Related  research  that  began  in  1996  resulted  in  the  NUVIS  instrument  of  today,  a 
hyperspectral  imager  with  a  spectral  bandwidth  of  300  to  375  nm  at  a  resolution  of 
approximately  0.3  nm.  (Cleary,  1998) 

A  hyperspectral  imager  produces  a  spectral  image  of  a  scene,  which  is  typically 
displayed  in  a  format  known  as  a  hyperspectral  cube.  Imaging  spectrometers  typically 
use  a  two-dimensional  array  (e.g.,  a  CCD)  and  produce  the  three-dimensional  data  cube 
containing  two  spatial  dimensions  and  a  third  spectral  dimension.  The  third  dimension  is 
obtained  by  scanning  either  the  spectral  dimension  or  one  of  the  two  spatial  dimensions. 
NUVIS  is  the  latter  type.  A  sample  of  NUVIS  output  is  demonstrated  in  the  image  cube 
illustrated  in  Figure  2.2.  A  large  number  of  images  are  created  from  contiguous,  rather 
than  disjoint  regions  of  the  spectrum.  This  technique  typically  yields  finer  resolution  than 
is  possible  with  a  multispectral  imager  and  in  doing  so,  provides  much  more  information 
to  the  data  analyst.  It  is  from  this  series  of  images,  each  measuring  the  level  of  incident 
energy  at  a  specific  wavelength,  that  a  combined  image  is  formed. 

Today,  countless  airborne  hyperspectral  sensors  are  employed  in  the  field.  Table 
2.2  provides  a  summary  of  many  of  these  instruments.  This  table  also  demonstrates  the 
relative  non-existence  of  UV  hyperspectral  imagers.  Thus,  NUVIS  is  a  portable,  ground 
based,  passive  remote  sensor  capable  of  hyperspectral  imaging  at  UV  wavelengths;  a 
prototype  instrument  that  could  someday  possibly  be  employed  as  an  airborne  sensor. 


9 


HYPERSPECTRAL 

CUBE 


Figure  2.  2.  NUVIS  hyperspectral  cube.  Image  illustrates  the  two  spatial  and  one  spectral 

dimension  comprising  the  image  cube. 
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Sensor 

(Agency) 

Number 
of  Bands 

Spectral 

Range(nm) 

Bandwidth 

FWHlVI(nm) 

IFOV 

(mrad) 

FOV 

(deg) 

Data 

Period  of 
Operation 

AAS  (ASTER) 

1 

760-850 

90 

1.0, 2.5, 

28.8,65 

Image 

Since(s) 

(DAIS-2815)  (GER) 

3 

3000-5000 

600-700 

or  5.0 

, 

Cube  (IC) 

1991 

20 

8000-12000 

200 

or  104 

AHS.Daedelus 

48 

440-12700 

20-1500 

2.5 

86 

IC 

1994 

AIS-2 

(NASA/JPL) 


AISA 

(Karelsilva  Oy) 


AMSS 

(Geoscan) 


Sensor 

(Agency) 


ARES  (AIP) 


ASAS 

Upgraded  ASAS 
(NASA/GSFC) 


AVIRIS  (JPL) 


CASI 

(Itres  Research) 


FTVHSI 

(Kestrel) 


GER-63 

Scanner 

(GER) 


1982-85 


Channel  24 
4 
29 
6 


400-1000 

1500-2000 

2000-2500 

8K-12.5K 


Continued  on  next  page. 
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Table  2.  2.  Summary  of  hyperspectral  airborne  imaging  spectrometers. 

From  Kramer,  1996. 


B.  ULTRAVIOLET  REMOTE  SENSING 
1.  The  Concept 

Passive  remote  sensing  at  ultraviolet  wavelengths  is  quite  dependent  on  the 
spectral  distribution  of  the  solar  irradiance  at  groimd  level.  This  distribution  is  illustrated 
in  Figure  2.3.  From  this  diagram,  it  is  easy  to  see  that  the  solar  irradiance  available  at 
ultraviolet  wavelengths  is  relatively  low  compared  to  that  at  visible  and  infrared 
wavelengths.  NUVIS  was  designed  to  operate  in  the  daytime,  using  ambient  scattered 
UV  radiation  in  the  300  to  400nm  spectral  band. 
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Figure  2.  3.  Spectral  distribution  of  the  solar  irradiance  (flux).  Notice  that  the  magnitude 
of  the  available  flux  falls  off  rapidly  below  400nm.  This  decline  is  due,  in  part,  to  ozone 
absorption,  which  will  be  discussed  later.  From  Finlay son-Pitts  &  Pitts,  1986. 


2.  Earth’s  Atmosphere  and  Ultraviolet  Wavelengths 

The  earth’s  atmosphere  can  be  divided  into  a  number  of  atmospheric  layers  each 
with  its  own  unique  characteristics.  Figure  2.4  illustrates  these  layers  and  provides  the 
approximate  behavior  of .  temperature  over  these  regions.  One  distinguishing 
characteristic  of  the  earth’s  atmosphere  is  that  it  is  generally  a  homogeneous  mixture  of 
nonreactive  major  component  gases  (N2,  O2,  noble  gases)  superimposed  upon  variable 
concentrations  of  minor  volatile  and/or  chemically  reactive  gases,  as  indicated  in 
Table  2.3. 
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Temperature  (K) 

Figure  2.  4.  Atmospheric  layers.  From  Lenoble,  1993. 


Constituent  Name  /  (type) 

Formula 

Molec.  Mass 

%  by  volume 

Nitrogen/  (Permanent) 

Nj 

28.0134 

78.084 

Oxygen  /  (Permanent) 

O2 

31.9988 

20.9476 

Argon  /  (Permanent) 

Ar 

39.948 

0.934 

Water  Vapor  /  (Variable) 

H2O 

18.0160 

0-7 

Carbon  Dioxide  /  (Variable) 

COj 

44.00995 

0.01-0.1 

Ozone  /  (Variable) 

03 

47.9982 

0-0.01 

Neon  /  (Permanent) 

Ne 

20.183 

0.001818 

Helium  /  (Permanent) 

He 

4.0026 

0.00052.4 

Methane  /  (Permanent) 

CH4 

16.04303 

0.0002 

Sulfur  Dioxide  /  (Variable) 

SO2 

64.064 

0-0.0001 

Hydrogen  /  (Permanent) 

H2 

2.01594 

0.00005 

Nitrogen  Dioxide  /  (Variable) 

issaBH 

46.0055 

0-0.000002 

Table  2.  3.  Atmospheric  constituents.  From  Durkee,  1999. 


Since  NUVIS  currently  operates  at  ground  level,  the  properties  of  ultraviolet 


radiation  in  the  lower  troposphere  are  of  paramount  interest.  Before  any  solar  ultraviolet 


radiation  can  reach  the  troposphere,  it  must  first  pass  through  the  various  overlying 
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atmospheric  layers.  Far  ultraviolet  wavelengths  are  strongly  absorbed  at  higher  altitudes 
(above  100km)  by  molecular  and  atomic  oxygen  and  nitrogen.  More  importantly,  a  great 
majority  of  UV  absorption  is  caused  by  oxygen  in  the  mesosphere  and  upper  stratosphere, 
principally  in  the  form  of  ozone.  These  absorbers  are  active  in  the  spectral  bandpass  of 
NUVIS  (-SOOnm).  Figure  2.5  summarizes  these  effects. 


Figure  2.  5.  Depth  of  penetration  of  solar  radiation.  Altitudes  correspond  to  an 
attenuation  of  1/e.  FromLenoble,  1993. 


C.  SULFUR  DIOXIDE 

1.  An  Atmospheric  Pollutant 

Absorption  of  ultraviolet  radiation  by  atmospheric  gases  is  not  limited  to  the 
cases  of  nitrogen,  oxygen  and  ozone  alone.  There  are  a  number  of  trace  gases  in  the 
atmosphere  that  absorb  at  these  shorter  (UV)  wavelengths.  Nitrogen  dioxide  (NO2) 
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and  sulfur  dioxide  (SO2)  are  two  such  examples.  For  the  purpose  of  this  work,  only  the 
trace  gas,  sulfur  dioxide,  will  be  studied. 

Sulfur  dioxide.  Figure  2.6,  is  a  colorless  gas,  which  has  a  characteristic 
suffocating  odor.  It  is  extremely  irritating  to  the  eyes  and  respiratory  tract.  It  is  formed 
in  the  combustion  of  sulfur-containing  fuels  and  in  the  smelting  of  sulfide  ores  and  is  a 
major  atmospheric  pollutant.  This  trace  gas  is  only  a  small  part  of  the  much  larger  sulfur 
cycle  illustrated  in  Figure  2.7.  This  cycle  not  only  demonstrates  that  SO2  is  a  primary 
pollutant  emitted  by  fuel  combustion  (from  coal  burning  power  plants)  but  also 
demostrates  that  most  volcanic  sulfur  is  emitted  as  SOj  (Finlayson-Pitts  and  Pitts,  1986). 


^Sulfur  (S) 

•  Oxygen  (O) 


Figure  2.  6.  Sulfur  dioxide  molecule  (SOj) . 


Figure  2.  7.  Sulfur  cycle  demonstrating  volcanic  and  man-made  sources  of  SO2. 
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Monitoring  of  sulfur  dioxide  has  become  a  global  concern  as  data  are  routinely 
collected  from  a  number  of  different  sensors  to  monitor  SO2  levels  world-wide. 
Monitoring  SO2  is  of  interest  for  a  number  of  reasons.  In  regions  of  high  volcanic 
activity,  SOj  monitoring  can  provide  yet  another  valuable  tool  to  assist  in  the  study  of 
volcanic  activity.  For  years  now,  the  Correlating  Spectrometer  (COSPEC),  a  non¬ 
imaging  spectrometer,  has  been  employed  by  vulcanologists  to  monitor  effluent  SO2 
levels  ( see  Figure  2.8  ).  Such  monitoring  is  important,  especially  in  densely  populated 
regions,  to  aid  in  predicting  volcanic  activity  and  to  warn  of  a  rise  in  SO2  concentrations 
at  ground  level. 


Figure  2.  8.  COSPEC  instrument  at  a  volcanic  site.  From  USGS,  1999. 
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The  Environmental  Protection  Agency  (EPA)  has  long  been  interested  in 
monitoring  anthropogenic  (man-made)  sources  of  sulfur  dioxide.  A  primary  focus  has 
been  placed  on  coal  burning  power  plants.  The  EPA  has  established  strict  guidelines  for 
allowable  types  and  levels  of  emissions  that  these  plants  can  produce.  Monitoring 
effluent  levels  from  these  sources  can  be  accomplished  either  by  sensors  installed  vdthin 
the  power  plant’s  stacks  themselves  {in  situ  measurements)  or  by  remote  sensing 
techniques. 

With  these  factors  in  mind,  NUVIS  designers  set  out  to  build  an  instrument 
capable  of  detecting  SOj  eniiissions  from  both  natural  and  man-made  effluent  sources. 
Since  typical  SO2  concentrations  range  from  only  tens  of  ppm  in  the  case  of  regulated 
power  plants  to  himdreds  or  even  thousands  of  ppm  for  volcanic  vents,  NUVIS  would  be 
required  to  detect  SO2  over  a  very  wide  range  of  absorption  column  densities. 

2.  An  Ultraviolet  Absorber 

Sulfur  dioxide  was  chosen  for  study  because  the  gas  exhibits  a  strong  absorption 
signature  in  the  near-ultraviolet  (200-400nm).  This  signature  is  illustrated  in  Figure  2.9. 
Although  SO2  absorbs  to  some  extent  over  the  entire  240  to  320nm  wavelength  range, 
the  region  longward  of  300nm  is  of  most  interest  for  measuring  SO2  concentrations  in  the 
atmosphere.  This  region  is  advantageous  because  the  illumination  source  for  doing 
absorption  spectroscopy  in  the  field  is  Rayleigh  scattered  solar  UV  radiation.  At  ground 
level,  this  scattered  radiation  diminishes  sharply  at  wavelengths  shorter  than  about  300nm 
(see  Figure  2.3).  A  further  advantage  is  that  if  in  this  region,  SO2  spectra  can  be 
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CROSS  SECTION  (cm^ 


readily  distinguished  from  the  spectra  of  other  absorbers  (see  Figure  2.10).  Hence,  we 
chose  to  concentrate  on  laboratory  and  field  measurements  of  SO2  spectra  in  the  300  to 
350  mn  wavelength  range. 


Figure  2.  9.  Gas  phase  absorption  spectrum  of  SO2.  The  extinction  coefficient  s  is 
proportional  to  the  absorptivity  coefficient  presented  in  section  I.D.2.  From  Atkins,  1990. 


Figure  2.  10.  Photoabsorption  cross  sections  of  atmospheric  UV  absorbers.  Solar  UV 
cross  section  data  obtained  from  the  Solar  Stellar  Irradiance  Comparison  Experiment 
(SOLSTICE)  is  also  displayed.  From  Strickland,  1999. 
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D.  ABSORPTION  SPECTROSCOPY 


1.  Introduction 

Absorption  spectroscopy  is  the  study  of  the  wavelength-dependent  absorption  of 
radiation  from  a  continuum  source  by  some  chemical  species  of  interest.  In  the 
laboratory,  this  technique  usually  involves  working  with  a  calibrated  continuum  light 
source.  In  NUVIS  field  applications,  the  background  continuum  source  is  the  ambient 
Rayleigh  scattered  UV  solar  radiation,  as  discussed  above. 

The  solar  irradiance  striking  the  atmosphere  has  approximately  the  spectral 
distribution  of  a  5800  K  blackbody  (see  Figure  2.3).  The  Raleigh  scattering  cross  section 
is  proportional  to  \IX^,  where  X  denotes  wavelength  (Lenoble,  1993).  The  combined 
effect  of  these  two  natural  phenomena  is  to  provide  a  useful  background  continuum 
source  in  the  300nm  wavelength  region.  The  COSPEC  instrument  previously  mentioned 
also  utilizes  this  continuum  radiation  source. 

Understanding  how  molecular  absorption  spectroscopy  is  used  to  measure 
molecular  concentrations  requires  a  basic  knowledge  of  the  energetics  of  the  molecules  of 
interest.  The  total  energy  £■  of  a  molecule  arises  from  four  contributions,  represented  by 
the  following  equation: 

E„ib  +  E„ec  +  E,„„5  .  (2- 1 ) 

E„t  is  the  rotational  energy  of  the  free  molecule  about  its  center  of  mass.  E^j^  is  the 
vibrational  energy  of  the  constituent  atoms,  typically  in  one  of  the  molecule’s  vibrational 
normal  modes.  E^i^^  is  the  electronic  energy  analogous  to  that  of  the  electronic  energy  of 
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isolated  atoms.  is  the  molecule’s  overall  translational  kinetic  energy.  In 

spectroscopic  applications,  the  rotational,  vibrational  and  electronic  energy  levels  are  of 
interest,  since  energy  is  quantized  for  these  modes  of  free  atoms  and  molecules.  Because 
of  the  bond  energies,  masses  and  moments  of  inertia  of  small  molecules,  transitions 
between  rotational  states  are  usually  observed  at  infrared  and  longer  wavelengths. 
Vibrational  transitions  are  typically  observed,  at  visible  and  infrared  wavelengths. 
Electronic  transitions  generally  occur  at  visible  and  UV  wavelengths.  (Lenoble,  1993) 

Figure  2.11  provides  a  generalized  representation  of  these  energy  levels  for  a 
simple  molecule.  Like  many  other  molecules,  gaseous  SO2  exhibits  its  own  unique 
characteristic  UV  absorption  spectrum.  If  an  incident  radiation  field  induces  photons 


Figure  2.  11.  Molecular  energy  curves  for  two  electronic  states.  These  curves 
demonstrate  the  relationship  between  vibrational,  rotational  and  electronic  levels. 
(This  general  depiction  is  not  to  scale.)  From  Eisberg  and  Resnick,  1985. 
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whose  energies  correspond  to  the  difference  between  two  energy  levels  of  the  molecule, 
the  molecule  may  absorb  a  photon  and  jump  to  a  higher  energy  level.  The  result  is  a 
wavelength-dependent  absorption,  or  absorption  spectrum,  that  is  characteristic  of  the 
quantum  structure  of  the  absorbing  species. 

2.  Terminology 

Transmittance  is  the  term  used  to  describe  the  fraction  of  the  energy  of  the 
incident  radiation  that  passes  through  an  absorbing  medium.  Figure  2.12  illustrates  a 
beam  of  parallel  radiation  before  and  after  it  has  passed  through  a  layer  of  absorbing 
medium  having  a  thickness  of  d  and  a  concentration  of  c  of  an  absorbing  species.  Due 
to  the  interactions  between  the  incident  photons  and  the  absorbing  particles,  an  incident 
intensity /(,  is  attenuated  to  an  exiting  intensity  /  (Skoog,  1985).  The  transmittance  T  of 


Absorbing  medium 
of  concentration:  ^ 


Path  length 


.  d 


Figure  2.  12.  A  parallel  beam  of  incident  radiation  before  and  after  it  has  passed  through 

a  layer  of  absorbing  .medium. 


the  medium  is  thus  defined  as  the  ratio 


Absorbance,  A  is  defined  as 


A=-log(I)=log(Io/I).  (2.3) 

If  the  likelihood  of  multiple  absorptions  and  reemissions  along  the  path  is  small,  there  is  a 
linear  relationship  between  absorbance  and  concentration  of  an  absorbing  species  given 
by  Beer’s  Law 

A=a(X)  d  c  (2.4) 

(Skoog,  1985).  The  fimction  a(X)  is  called  the  absorptivity  coefficient.  From  Beer’s 
Law,  the  attenuation  of  a  signal  can  be  given  by 

/=/(,  exp(-(jN  d),  (2.5) 

where  cris  the  absorption  cross  section  and  iVis  the  molecular  number  density. 

Absorption  is  a  quantity  that  differs  from  absorbance:  It  can  be  described  in  terms 
of  the  Einstein  probability  coefficient,  B12,  but  in  practice,  it  is  more  often  described  in 
terms  of  an  absorption  coefficient,  k(v).  This  parameter  describes  the  fractional  decrease 
in  flux  density  (intensity)  at  frequency  v  per  unit  path  length  through  the  absorbing 
medium: 

-5I(v)=J(v)  k(v)  5k,  (2.6) 

where  I( v)  is  incident  upon  the  layer  of  thickness  &c.  Assuming  a  homogeneous  layer  of 
thickness  d,  integration  yields 

J/^)  =  Io( y)  exp(-k( v)  d),  (2.7) 
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where  k(v)  d  is  known  as  the  frequency-dependent  optical  depth,  nf v).  The  effect  of 
increasing  optical  depth  on  an  isolated  absorption  line  profile  is  illustrated  in  Figure  2.13. 
(Thome  et  al.,  1999) 

From  a  practical  standpoint,  it  is  the  effective  area  of  an  absorption  line  that  gives 
the  most  useful  measure  of  spectral  absorption  that  is  taking  place.  The  areas  under  the 
absorption  lines  (curves)  of  Figure  2.13,  may  be  defined  by  the  following  equation. 


Figure  2.  13.  Hypothetical  absorption  line  profiles  demonstrating  the  effect  that  an 
increase  in  optical  depth  has  on  a  line  shape.  Lines  a,  b,  and  c  represent  shapes  for 
increasing  optical  depths,  T(y)  .  From  Thome  et.  al.,  1999. 


which  is  normalized  to  the  incident  flux  density. 


W=fu„e  (l-expf-k(v)J)  dv.  (2.8) 
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Equation  2.8  equates  well  to  transmittance  through  the  following  relation, 

W=Ul-T(v)Idv,  (2.9) 

where  W  is  called  the  line’s  equivalent  width,  and  T(v)  is  the  frequency  dependent 
transmittance  from  equation  2.2.  Figure  2.14  further  depicts  the  equivalent  width  of 
spectral  line  where  the  shaded  area  is  equivalent  to  the  area  between  the  line  profile  and 
the  continuum  baseline,  1(0). 


Figure  2.  14.  Equivalent  -width  measurement.  Note  that  the  shaded  area  is  equal  to  the 
area  between  the  line  profile  and  continuum  baseline,  1(0).  From  Thome  et.  al.,  1999. 

It  should  be  noted  that  the  actual  absorption  of  an  isolated  spectral  line  can  be 
spread  over  a  broad  range  of  frequencies  by  natural,  Doppler  and  pressure  broadening 
(Thome  et  al.,  1999).  Furthermore,  the  relation  between  the  absorbance  and  I(v)  is 
logarithmic,  not  linear,  so  equation  2.9  does  not  hold  tme  for  all  cases.  Measuring 
individual  equivalent  line  widths  is  thus  a  more  detailed  process  than  measuring  the 
overall  trend  of  the  absorption  cross  section  of  Equation  2.5.  However,  it  can  be  shown 
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that  linearity  applies  as  long  as  the  absorbing  line  is  not  highly  saturated  ( »  0  at  the 
profile’s  center  ).  Therefore,  for  the  purposes  of  this  thesis  research.  Equation  2.9  was 
utilized  in  calculating  equivalent  widths  using  nonsaturated  lines.  This  was 
accomplished  without  correcting  for  any  line  shape  effects  and  assumed  that  Equation 
2.9  was  applicable  to  all  equivalent  width  measurements.  For  a  more  detailed  approach, 
see  Thome  et.  al,  1999. 

3.  The  Curve  of  Growth 

In  spectrochemical  analysis  by  means  of  absorption  spectroscopy,  it  has  become 
standard  practice  to  quantify  the  effects  of  an  absorbing  species  abundance  on  the  profile 
of  an  absorption  line  as  an  empirical  graph  called  a  ciuve  of  growth.  The  curve  of 
growth  for  a  spectral  line  describes  the  behavior  of  the  equivalent  width  as  a  function  of 
concentration  of  the  absorbing  species.  It  normally  takes  the  form  of  a  log-log  plot  as 
illustrated  in  Figure  2.15.  In  practice,  the  equivalent  widths  are  measured  for  a  series  of 
known  absorber  concentrations  in  the  laboratory.  The  equivalent  width  of  lines  found  in 


Figure  2.  15.  Typical  curve  of  growth.  (Also  known  as  a  calibration 

curve.) 
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field  data  may  be  subsequently  compared  to  the  curve  of  growth,  which  then  yields  the 
molecular  abundance  of  the  field  sample. 

Remarkably,  curve  of  growth  analysis  has  been  found  to  be  relatively  insensitive 
to  different  line  profile  shapes  (Thome  et.  al.,  1999).  Hence,  it  is  of  general  utility  across 
a  wide  variety  of  spectroscopic  instruments  and  conditions.  Two  commonly  encountered 
line  shapes,  Gaussian  and  Lorentzian,  are  illustrated  in  Figure  2. 16.  This  study  will  only 


Figure  2.  16.  Gaussian  and  Lorentzian  profiles.  Areas  under  each  curve  are  the  same. 

From  Thom  etal.,,  1999. 

involve  calculations  based  on  Gaussian  line  profiles  that  may  be  represented  by  curves 
having  the  generalized  form: 

F(x)=Ao  exp{-f(x-Aj)/A2f/2},  (2.10) 
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where  F(x)  is  the  profile,  Agxs  the  height  of  the  Gaussian,  Aj  is  the  center  of  the  Gaussian 
and  ^2  is  the  width  of  the  Gaussian.  Equation  2. 1 0  will  be  utilized  in  the  curve  of  growth 
calibration  of  the  NUVIS  instrument  and  in  the  analysis  of  field  data  presented  in 
Chapters  IV  and  V. 
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III.  NUVIS 


A.  BACKGROUND 

NUVIS  was  designed,  developed  and  constructed  by  former  NPS  Associate 
Professor  Dr.  David  Cleary  and  students,  Todd  Hooks  (Hooks,  1997)  and  Andrew 
MacMannis  (MacMannis,  1997).  It  is  the  latest  in  a  succession  of  ultraviolet 
spectroscopic  instruments  developed  by  Dr.  Cleary  and  is  the  first  to  incorporate  spectral 
imaging  capabilities  (Cleary,  1998).  For  a  more  detailed  account  of  some  of  this  earlier 
work  done  on  UV  instruments  such  as  MUSTANG  and  DUUVIS,  see  Johnson,  1996. 

B.  OPTICAL,  OPTOMECHANICAL  AND  ELECTRONIC  HARDWARE 

In  principle,  the  NUVIS  instrument  of  today  is  identical  to  the  instrument  first 

developed  by  Hooks  and  MacMannis.  Designed  and  built  to  test  the  utility  of 
hyperspectral  imaging  at  ultraviolet  wavelengths,  the  instrument’s  optical  design  is  a 
classic  Rowland  circle  spectrograph  and  its  imaging  design  follows  that  of  a  push  broom 
scaimer.  As  it  appears  in  Figure  3.1,  NUVIS  is  a  portable  instrument  capable  of  field 
operation  with  ease.  Instrumental  parameters  remain  as  originally  detailed  by  Dr.  Cleary 
and  are  summarized  in  Table  3.1. 

The  instrument  itself  is  comprised  of  nine  basic  components,  which  are  illustrated 
in  Figure  3.2.  Light  enters  the  instrument  through  a  filter  window,  designed  to  pass  only 
UV  wavelengths.  It  then  encounters  a  scanning  mirror  assembly.  This  assembly  is 
comprised  of  three  sub-components:  1)  the  mirror  itself,  2)  a  stepper  motor  driver,  which 
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Figure  3.1.  NUVIS  instrument. 


NUVIS  PARAMETERS 

Optical  Design 

Rowland  Circle  Spectrograph 

Imaging  Design 

Push  Broom  Scanner 

Spectral  Response 

300  to  375nm 

Spectral  Resolution 

0.3  nm 

Full  Field  of  View 

20'’x6“  (350mradx  lOOmrad) 

Angular  Resolution 

~  Imrad 

Bands  (Spexels) 

640 

Positions 

480 

Lines  (available) 

351>ZZZ 

Mode  of  Operation/ 

SNAP/  640x480x1  data  array 

DATA  output 

SCAN/  image  cube  (640  x  480  x  ZZZ) 

Table  3.1.  NUVIS  parameters.  From  Cleary,  1998. 
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Figure  3.2.  NUVIS  components.  Top:  Schematic  of  NUVIS:  1)  filter  window,  2) 
scanning  mirror  assembly,  3)  baffles,  4)  telescope  mirror,  5)  slit,  6)  diffraction 
grating,  7)  image  intensifier,  8)  CCD  camera,  and  9)  cable  coimectors. 
Bottom:  Photo  of  NUVIS  instrument  components,  (with  cover  removed) 
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rotates  the  mirror,  and  3)  the  shaft  encoder,  which  provides  feedback  information  about 
the  angular  position  of  the  mirror.  Light  reflected  by  the  scanning  mirror  strikes  a 
telescope  mirror,  which  focuses  the  beam  onto  the  entrance  slit,  which  is  located  on  the 
Rowland  circle.  From  the  slit,  light  travels  to  a  field-flattened  spherical  diffraction 
grating.  The  diffracted  light  then  goes  to  an  image  intensifier.  The  image  intensifier, 
optimized  for  UV  operation,  optically  amplifies  the  signal  by  a  factor  of  nearly  1 0®.  The 
image  intensifier’s  output  is  coupled  to  a  Pulnix  TM-745e  charge  coupled  device  (CCD) 
video  camera.  The  camera  has  768  by  494  pixel  resolution  with  1  Ipm  by  ISpm  pixels. 
The  camera  output  is  conventional  NTSC  30  frames  per  second  (tps)  analog  video 
through  a  standard  RS-170  interface  and  is  limited  to  an  eight  bit  grayscale  and  a 
resolution  of  640  by  480  lines.  (Hooks,  1997) 

C.  THE  COMPUTER,  DATA  ACQUISITION  AND  CONTROL  SYSTEMS 

The  next  major  component  of  the  total  NUVIS  system  is  the  NUVIS  data 
acquisition  and  control  computer,  including  the  CPU,  monitor  and  supporting  peripherals. 
The  233  MHz  PENTIUM  Pro  PC  includes  256Mb  Ram  and  a  3.0  GB  hard  drive.  It  is 
equipped  with  standard  3.5  floppy  and  CD-ROM  drives  as  well  as  a  1  GB  internal  JAZ 
drive.  The  video  output  from  NUVIS’s  camera  is  processed  into  digital  form  by  a  PC 
IMAQ-1408  frame  grabber  interface  card.  Interfaces  to  the  stepping  motor  controller  and 
shaft  position  encoder  are  provided  by  DAC-02  and  PIO-24IS A/EISA  cards.  The 
monitor  is  an  NEC  Multisync  LCD  15  inch  flat  screen  display.  Figure  3.3  illustrates 
external  views  of  the  computer  system. 
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Cl'U,  Monitor,  Keyboard,  ana  mouse.  Note  jtsacK  panel  oi  uru  witn  networK  nub 
that  all  components  are  easily  operated  while  attached.  DAC-02  ,  PIO-24ISA/EISA  &  the 
still  inside  the  customized  shipping  case.  PC  IMAQ-1408  frame  grabber  cards 

displayed  from  the  bottom  up- 

Figure  3.3.  NUVIS  control  computer. 

Control  of  the  scanning  mirror  drive  mirror  drive  is  provided  by  the  MAX-410 
stepper  motor  controller  that  is  illustrated  in  Figure  3.4.  Figure  3.5  shows  the  complete 
NUVIS  system  deployed  in  the  field,  demonstrating  some  of  the  logistics  and  setups 
involved  with  on-site  operations.  A  detailed  account  of  suggested  system  components 
required  for  field  operation  can  be  found  in  the  NUVIS  field  checklist  of  Appendix  A. 


Figure  3.4.  Max-410  stepper  controller. 
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Left:  Volcanic  Site 


Right:  Power  Plant  Site 


Figtore  3.5.  NUVIS  system  components  and  field  logistics, 
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D.  DATA  ACQUISITION,  CONTROL  AND  DATA  DISPLAY  SOFTWARE 

As  NUVIS  was  being  develpoed,  several  key  software  routines  were  developed 
for  instrument  control,  data  acquisition  and  data  display.  Various  modules  of  this 
software  were  written  in  Microsoft  Visual  Basic  and  Interactive  Data  Language  (IDL), 
which  is  a  Research  Systems,  Inc.  (RSI)  product.  Specifics  of  the  software  may  be  found 
in  the  following  references:  Hooks,  1997;  Aitken,  1997;  and  Using  IDL,  1997  as  well  as 
in  Appendix  B. 

In  essence,  the  acquisition  of  spectral  image  data  with  NUVIS  is  a  simple 
algorithmic  loop.  The  scanning  mirror  is  rotated  so  that  one  edge  of  the  target  scene  of 
interest  falls  on  the  instrument’s  entrance  slit.  The  diffraction  grating  then  disperses  the 
light  from  the  slit  across  the  camera.  Because  the  Rowland  circle  spectrometer  is 
stigmatic  (image-preserving),  the  resultant  2D  image  consists  of  spectral  information 
along  the  dispersion  axis  as  a  function  of  ID  scene  height  along  the  orthogonal  axis.  In 
other  words,  an  image  frame  consists  of  a  set  of  spectra  for  each  scene  pixel  focused  on 
the  entrance  slit.  Each  of  these  scene  pixels,  however,  belongs  to  only  a  narrow  vertical 
strip  of  the  whole  target  scene.  So,  the  camera  frame  is  stored,  then  the  scanning  mirror 
is  rotated  slightly  to  position  the  next  vertical  strip  of  the  scene  onto  the  exit  slit,  and  the 
frame  acquisition  procedure  is  repeated.  This  sequence  of  operation  is  repeated  until 
data  from  the  entire  scene  has  been  acquired.  Subsequent  off-line  processing  can  then  be 
used  to  generate  a  hyperspectral  cube  data  set  and  to  display  the  data  in  various  formats. 
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Sample  output  from  a  field  deployment  is  provided  in  Figures  3.6  through  3.8. 
This  particular  set  is  drawn  from  data  taken  of  an  industrial  coal  burning  power  plant. 
The  image  scene  includes  three  operating  smokestacks.  All  three  stacks  are  producing 
some  level  of  sulfur  dioxide  emissions,  as  well  as  a  number  of  other  effluent  products. 
Two  of  the  stacks  have  modem  to  contain  SO2  scmbbers  that  reduce  emission  levels.  The 
third  (dirty)  stack  emits  SO2  at  levels  on  the  order  of  ten  times  the  ppm  values  of  the 
scmbbed  stacks.  This  particular  collection  was  conducted  on  site  at  a  range  of  1.5  km. 
(Further  detail  of  this  field  collection  may  be  found  in  Chapter  V.) 

The  black  and  white  image  of  Figure  3.6  is  produced  directly  from  the  sensor’s 
recorded  digital  signals  through  a  transpose,  rotation  and  resizing  of  the  data  array. 
Although  the  stacks  are  quite  visible,  this  format  does  not  provide  enough  information  to 
allow  for  the  isolation  and  identification  of  SO2  signatures  in  the  data.  A  false-color 


Figure  3.6.  NUVIS  black  and  white  image.  This  image  was  produced  by  the 
IDL  procedure,  Bw_roughimage.pro.  From  Cleary,1998. 
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image  can  also  be  created  to  survey  the  data.  This  was  done  for  the  power  plant  work, 
using  color  assignments: 


Red  :  350-375nm, 

Green:  325-350ntn, 

Blue :  300-325nm. 

This  approach  results  in  the  false  color  image  of  Figure  3.7,  in  which  potential  regions  of 
sulfur  dioxide  appear  yellow  in  color. 


Figure  3.7.  NUVIS  false-color  image.  This  is  a  product  of  the  IDL  procedure, 
"color_image.pro".  Regions  of  SO2  appear  yellow  in  color.  From  Cleary,  1998. 


Field  operation  of  NUVIS  necessitates  a  means  of  obtaining  quick  looks  at  the 
data  in  order  to  ascertain  whether  or  not  the  instrument  is  functioning  properly.  An  IDL 
procedure  was  developed  to  allow  thus.  It  performs  a  least-squares  fitting  of  the  data  to  a 
modeled  absorption  cross  section  of  SO2.  Output  displays  from  this  procedure  are 
illustrated  in  Figure  3.8.  This  corresponds  to  200  different  spatial  (altitude)  vs.  spectral 
(wavelength)  slices  of  the  hyperspectral  image  cube.  For  presentation  purposes,  data 
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processing  has  been  interrupted  midway  through  the  scene  at  line  number  99  to  illustrate 
the  output  windows  presented  to  the  data  analyst.  Window  1  shows  a  display  similar  to 
that  of  Figures  3.6  and  3.7.  During  the  execution  of  this  procedure,  processing  is 
actively  displayed  from  from  right  to  left.  A  full  false-color  image  of  the  scene  is  first 
presented  and  is  then  actively  processed,  scan  by  scan  (line  by  line)  within  the  window. 
The  display  labeled  as  Window  2  represents  the  raw  data  or  digital  number  that  is 
recorded  for  each  individual  pixel  in  the  spatial  (altitude)  vs.  spectral  (wavelength)  slice 
(plane),  corresponding  to  line  number  99.  The  raw  data  displayed  on  the  right, 
corresponds  to  the  line  dividing  the  color  and  black  and  white  images  in  the  window  on 
the  left.  This  line  marks  the  point  at  which  the  image  cube  was  sliced  to  produce  such  a 
display.  Within  Window  2,  the  darker  pixels  represent  lower  levels  of  detection.  The 
opposite  is  true  for  the  brighter  pixels. 

Window  2  also  parallels  the  actual  raw  data  display  that  is  presented  to  the  user 
during  instrument  operation.  This  raw  display  is  the  only  real-time  feedback  provided 
during  system  operation.  Feedback  is  provided  through  the  NUVIS  operational  window, 
which  is  displayed  by  the  Visual  Basic  instrument  control  program  NUVIS.EXE.  A  view 
of  the  same  data  at  a  time  later  in  analysis  process  is  presented  in  the  lower  windows  one 
and  two.  Thus,  Figure  3.8  also  illustrates  the  representative  output  at  the  150*  vertical 
scan  line. 

In  the  field,  the  data  acquisition  and  control  PC  can  be  networked  to  another 
computer,  typically  a  laptop  PC.  The  various  IDL  data  analysis  and  display  procedures 
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can  then  be  run  on  the  laptop  while  the  host  PC  is  otherwise  occupied  with  running  the 
instrument.  This  provides  a  quasi-real-time  means  of  verifying  NUVIS’s  operation  in  the 
field. 


Window  1  X - ►  Window  2  A,  - ► 


Window  1  X - ►  Window  2  A  - ► 


Figure  3.8.  SO2  algorithm  output.  Upper  Windows  demonstrate  output  displays 
for  the  99*  vertical  line  scan  of  the  image  scene  during  execution  of  the  algorithm, 
"assay_image.pro."  Lower  windows  demonstrate  processing  some  time  later  at 
the  150*  scan  position.  In  each  case,  altitude  is  displayed  in  the  vertical  (spatial- 

Y)  direction.  From  Cleary,  1998. 
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IV.  LABORATORY  CALIBRATIONS  OF  NUVIS 

A.  PURPOSE 

In  order  for  NUVIS  to  be  applied  to  field  measurements  of  SO2  plumes,  two  types 
of  calibrations  must  be  performed,  either  beforehand  or  on  site  during  field  deployment. 
First,  an  approximate  spectral  wavelength  scale  must  be  established,  so  that  spectral 
signatures  of  interest  can  be  reliably  and  reproducibly  identified.  Second,  the  relevant 
curves  of  growth  for  the  signatures  must  be  available  so  that  the  instrument  can  be 
employed  for  reliable  and  credible  spectrochemical  abundance  measurements.  This 
chapter  summarizes  the  laboratory  calibrations  that  were  performed  in  these  contexts. 
The  next  chapter  discusses  applications  of  the  calibrations  to  actual  field  data. 

B.  WAVELENGTH  CALIBRATION 

Before  continuing  onto  generation  of  the  curves  of  growth,  it  is  necessary  to  first 
discuss  the  association  of  instrument  spectral  samples  (sometimes  known  as  spexels)  with 
wavelength.  This  phase  of  calibration  was  conducted  by  Dr.  Cleary,  just  prior  to  my 
involvement  with  the  NUVIS  project.  Results  of  the  wavelength  calibration  were  saved 
in  a  1 .2  MB  data  file  entitled  wl_cal.dat. 

The  calibration  was  accomplished  using  a  standard  spectroscopic  emission  line 
source.  The  source  used  was  a  platinum  (Pt)  hollow  cathode  lamp,  whose  characterisitics 
and  calibration  are  traceable  to  the  National  Institute  of  Standards  and  Technology 
(NIST).  This  source  was  chosen  because  it  exhibited  a  number  of  useful  spectral 
emissions  in  the  300  to  400nm  wavelength  range.  This  particular  lamp  is  used  in 
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conjunction  with  a  Harrison  6522  DC  power  supply.  For  operation,  the  supply  is  initially 
set  at  400  volts  and  0mA.  The  Current  is  gradually  increased  to  the  operational  level  of 
1 5mA  Nine  Pt  emission  lines  were  identified  and  utilized  in  the  calibration.  These  lines 
are  listed  in  Table  4. 1 . 


Line 

Wavelength 

Line 

Wavelength 

Line 

Wavelength 

1 

304.26 

4 

337.82 

7 

359.36 

2 

306.47 

5 

340.81 

8 

369.42 

3 

334.54 

6 

352.05 

9 

372.71 

Table  4. 1 .  Nine  platinum  emission  lines  chosen  for  the  wavelength  calibration. 

From  Cleary,  1998. 


As  with  any  diffraction  grating  instrument,  the  grating  equation  applies.  (For 
further  explanation  see  both  MacMannis,  1997  and  Hecht,  1998.)  This  means  that  the 
actual  wavelength  of  a  spectral  feature  should  be  a  sinusoidal  function  of  its  position 
(spectral  sample)  in  the  measured  data  frame.  It  was  found  experimentally,  however,  that 
the  sinusoid  could  be  adequately  approximated  by  a  linear  fit  across  the  300  to  400  nm 
range.  This  behavior  is  summarized  in  Equation  4.1 

Y=(0.118)  X  +  300.775,  (4.1) 

where  x  is  the  spexel  number,  m  the  spectral  dispersion,  b  the  wavelength  of  the  “zeroth” 
spexel  and  Y  is  the  corresponding  wavelength  in  nanometers.  This  behavior  is  also 
illustrated  in  Figure  4.1.  This  calibration  was  applied  in  all  instances  where  it  became 
necessary  to  associate  wavelength  to  a  corresponding  spectral  sample.  It  also  provided  a 
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simple  means  for  calibrating  the  wavelength  scales  of  spectra  measured  with  NUVIS, 
provided  that  the  optomechanical  configuration  of  the  instrument  remained  stable. 


Wavelength  Calibration 


Figure  4.1 .  Wavelength  calibration  associating  band  number/spexel  to 

wavelength. 


C.  EXPERIMENTAL  OVERVIEW 

As  stated  earlier,  only  absorption  lines  existing  longwards  of  305nm  were  chosen 
for  laboratory  measurement  because  absorption  lines  that  appeared  in  the  field  data  at 
shorter  wavelengths  were  either  saturated  or  subject  to  severe  noise  because  of  rapidly 
decreasing  background  illumination.  Figure  4.2,  vividly  demonstrates  the  relatively 
closely  spaced  group  of  absorption  lines  that  were  targeted  for  calibration.  In  the  end, 
four  absorption  lines  were  identified  as  suitable  for  equivalent  width  measurements  in  the 
calibration.  The  spectral  sample  numbers  and  wavelengths  corresponding  to  the  centers 
of  each  of  these  absorption  lines  are  listed  in  Table  4.2.  Because  NUVIS  was  expected  to 
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operate  over  a  wide  range  of  column  densities,  the  decision  was  made  to  first  concentrate 


on  the  instrument’s  lower  sensitivity  limits. 


Measured  S02  Absorption  Lines 


0  100  200  300  400  500  600  700 

Spectral  Sample  Number 

Figure  4.2.  Example  of  S02’s  closely  spaced  absorption  lines.  In  this 
plot,  lower  spectral  sample  numbers  equate  to  higher  UV  wavelengths 
and  higher  numbers  equate  to  lower  ones.  See  Figure  4. 1 . 


Measured 
Absorption  Line 
# 

1 

2 

3 

4 

Spectral  Sample 
Number 

535 

554 

572 

590 

Center 

Wavelength 

(nm) 

313 

310.9 

308.8 

306.5 

Table  4.2.  Absorption  lines  targeted  for  measurement. 


D.  LABORATORY  MEASUREMENTS  OF  SO2 

There  are  two  contraints  on  the  gas  concentrations  that  are  useful  for  laboratory 
calibration  studies.  First,  the  experimental  noise  floor  sets  a  practical  lower  limit  to  the 
instrument’s  sensitivity.  That  is,  weak  spectral  features  that  are  buried  in  the  noise  must 
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simply  go  unmeasured.  This  translates  into  a  lower  useful  limit  to  SO2  concentrations 
used  in  the  laboratory.  Second,  as  discussed  earlier,  highly  saturated  spectral  features  do 
not  provide  reliable  measurements  of  spectrochemical  abundance,  and  so  should  be 
avoided.  This  sets  a  practical  upper  limit  on  laboratory  gas  concentrations.  After 
weighing  such  considerations,  six  SO2  column  abundances  were' decided  upon  for  study. 
They  are  summarized  in  Table  4.3. 


Sample 

Column  abimdance 
(cc  concentration) 
[(#molecules/m^)  m)] 

lA 

5.33  X  10'* 

2B 

5.61  X  10'* 

3A 

7.48x10'* 

3A 

1.07  X  10" 

4B 

2.25  X  10" 

4B 

4.24x10" 

Table  4.3.  The  six  test  cell  concentrations  utilized  in  generating  the  curves  of 

growth. 


1.  Gas  Test  Cells 

SO2  was  placed  in  transparent  test  cells  of  various  lengths  and  illuminated  with 
UV  continuum  radiation  to  produce  laboratory  spectra.  The  first  test  cell  employed  was  a 
1.55  meter  long  PVC  tube  with  quartz  windows  (see  Figure  4.3).  Once  evacuated  with  a 
conventional  fore  pump,  the  cylinder  could  be  filled  through  a  pressure  regulator  with 
0.11  %  sulfur  dioxide  in  an  inert  dilution  gas  (N2).  The  valve  manifold  and  Matheson 
pressure  gage  ( 0  to  760  mm  Hg  )  that  were  both  utilized  in  this  process  are  illustrated  in 
Figure  4.4.  For  various  reasons,  the  1.55  m  cell  proved  to  be  cumbersome,  prone  to 
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leaks,  and  too  long  for  producing  low  column  abimdances  of  SOj.  Therefore,  another 
test  cell  was  designed. 


Figure  4.4.  1.55m  gas  test  cell  #1 . 
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Figure  4.5.  Gas  test  cell  #2. 

To  remedy  these  differences  with  the  first  two  gas  test  cells,  it  was  decided  to 
obtain  a  more  enriched  mixture  of  sulfur  dioxide  (Matheson  10%  mixture  in  air)  and  to 
use  much  shorter  gas  test  cells.  Four  inch  diameter  (10.16cm)  by  3.25cm  inner  path 
length  quartz  cells  were  purchased  from  Weiss  Scientific  Glass  Blowing  Company.  They 


The  next  test  cell  was  considerably  shorter  and  more  portable,  and  could 
potentially  be  utilized  not  only  in  the  laboratory  but  also  in  the  field.  It  is  illustrated  in 
Figure  4.5.  Even  though  this  design  was  more  portable  than  the  first,  it  was  quickly 
abandoned  for  two  reasons:  1)  The  test  cell  could  not  reliably  be  coupled  to  the  optical 
entrance  (filter  window)  of  NUVIS.  2)  The  cell  combination  was  mechanically  unstable 
and  could  potentially  damage  the  instrument’s  filter  window. 
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are  shown  in  Figure  4.6.  These  smaller  cells  also  had  a  mechanical  advantage:  they  can 
be  mounted  easily  at  NUVIS’s  optical  entrance  filter  window,  permitting  easy  laboratory 
and  field  use.  A  minor  difficulty  with  their  use  was  discovered,  however.  Their  available 
window  diameter  is  only  8.77cm  once  the  cell  is  installed  in  the  shop  fabricated 
aluminum  housing/mounting.  This  cause  slight  vignetting  at  the  comers  of  the 
instrument’s  field  of  view.  During  thesis  research,  care  was  taken  to  exclude  these 
regions  and  a  new  set  of  six  inch  diameter  test  cells  and  associated  mountings  are  under 
development.  In  Figure  4.7  a  layout  of  the  filter  window  and  the  gas  test  cells  is  shown. 


Figure  4.6.  Four  inch  gas  test  cells.  Photos  taken  before  and  alter  installation. 


-03.450 


Figure  4.7.  Test  cell  comparison.  Dotted  line  represents  NUVIS  filter 
window.  Circles  represent  available  cross  sections  associated  with  the  4  inch 
(inner  circle)  and  6  inch  (outer  circle)  cells  respectively.  The  reduced  cross 
sections  are  due  to  the  cell  mounts  themselves. 
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The  filter  window  is  represented  by  a  dotted  line  and  the  two  circles  (of  diameter  3.45 
and  5.475  inches,  respectively)  represent  the  available  cross  sections  provided  by  the  four 
inch  and  new  six  inch  cells. 

2.  Calibration  Source  and  Reflective  Screen 

During  initial  calibration  attempts,  an  EG&G  Gamma  Scientific  deuturium  lamp 
was  utilized  as  a  background  continuum  illumination  source.  In  the  300  to  375  nm 
spectral  region,  this  lamp’s  maximum  specific  output  irradiance  was  on  the  order  of 
10'^  (pW/cm^  nm).'  These  source  levels  proved  to  be  too  low  for  practical  absorption 
measurements  in  the  laboratory.  A  1000  Watt  quartz  halogen  tungsten  filament  lamp 
(FEL-309  IR)  from  Optronics  Laboratories  was  chosen  to  replace  the  deuterium  lamp.  It 
is  a  NIST  tracable  UV  source  and  was  utilized  throughout  all  calibrations.  The  spectrum 
of  this  source  is  illustrated  in  Figure  4.8. 


Spectral  Irradiance  of  F-309  at  50cm 


Figure  4.8.  Spectral  irradiance  of  FEL  calibration  source  over  operational 

NUVIS  wavelengths. 
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The  lamp  was  mounted  in  an  enclosed  box  with  a  controllable  output  aperture, 
which  was  used  to  maximize  the  optical  fluence  on  a  diffuse  reflective  screen,  described 
below.  The  lamp  was  powered  by  a  Hewlett-Packard  6030-A  power  supply.  Maximum 
drive  ciurent  was  set  to  8.0  amps.  The  lamp  required  a  15  minute  warm  up  period 
(according  to  the  manufacturer)  to  achieve  steady  state  operation  before  measurements 
could  be  taken. 

The  light  from  the  PEL  lamp  wps  scattered  by  a  reflective  screen  assumed  to  be 
an  ideal  diffuse  Lambertian  surface  (see  Hecht,  1998).  The  screen  was  a  Spectralon 
diffuser  with  greater  than  97%  reflectivity  in  the  wavelength  region  of  interest 
(  see  Figure  4.9  ).  It  was  placed  150  cm  from  the  PEL  and  approximately  normal  to  the 


Spectralon  Reflectance  Target  Calibration 


Figure  4.9.  Screen  reflectance  values. 
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lamp-screen  axis.  The  NUVIS  instrument  was  oriented  at  45  degrees  from  this  normal 
axis.  The  geometry  of  the  setup  is  illustrated  in  Figure  4.10. 


Figure  4. 1 0.  Experimental  set  up.  Instrument  is  placed  at  45  degree  angle 
from  screen-lamp  axis.  UV  source  is  enclosed  by  the  box  and  is 
positioned  for  normal  incidence  to  the  screen. 


3.  Test  Cell  Preparation 

Calibration  cells  were  prepared  by  first  connecting  them  to  the  valve  manifold.  A 
vacuum  was  drawn  on  the  cell  using  a  mechanical  fore  pump  vented  through  a  gas 
exhaust  hood.  The  cell  gas  pressure  was  monitored  by  the  Matheson  pressure  gage.  For 
background  calibrations,  evacuated  cells  could  be  used.  If  the  cell  was  to  be  filled  with 
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sulfur  dioxide,  the  appropriate  concentration  was  introduced  by  using  a  standard  gas 
regulator  while  monitoring  the  pressure. 

The  entire  setup  (see  Figure  4.11)  had  some  inherent  limitations.  The  gage  only 
marked  pressures  in  increments  of  5mm  of  mercury  (mm  Hg).  This  did  not  permit  very 
precise  determination  of  cell  gas  pressure  and  caused  uncertainties  in  gas  concentration, 
particularly  at  low  pressure. 


Figure  4.1 1.  Cell  preparation  setup. 


4.  Data  Acqusition 

Once  a  cell  was  prepared,  it  was  mounted  in  front  of  the  filter  window  of  the 
instrument.  Because  flourescent  lights  exhibit  well  known  emissions  in  the  ultraviolet, 
the  lab  was  darkened  before  measurements  were  taken.  A  range  of  instrument  gains  was 
first  tested  to  determine  where  sensor  saturation  occurred.  This  was  accomplished  by 
using  the  VIDEO  feature  displayed  in  the  NUVIS  control  window.  An  image  intensifier 


52 


gain  of  7.0  was  set  in  the  control  window.  Measurements  were  obtained  by  first  setting 
the  number  of  images  per  SNAP  to  100.  Spectrum  were  then  recorded  using  the  NliVIS 
SNAP  mode,  in  which  the  scan  mirror  is  not  rotated.  In  addition  to  gas  spectra  at  each 
concentration,  spectra  of  empty  cells  were  recorded  for  comparison  purposes.  The 
spectra  were  saved  to  a  1  GB  Jaz  cartridge  as  a  30MB  NUVIS  file. 

5.  First  Results 

For  the  first  experiment,  fifteen  different  spectra  were  recorded.  Three  representative 
spectra  from  experiment  #1  are  illustrated  in  Figure  4.12.  These  plots  represent  mean 
spectra  calculated  by  averaging  over  all  480  pixels  in  the  vertical  image  column  and  then 
averaging  over  all  1 00  images  to  enhance  the  experimental  ratio.  Only  spectral  samples  > 
500  are  displayed.  As  discussed  earlier,  absorption  is  characterized  by  the  ratio  ///<,. 
Hence  each  SOj  spectrum  was  then  divided  by  the  corresponding  spectrum  of  the 
evacuated  cell.  The  results  are  illustrated  in  Figure  4.13.  Neither  of  these  initial 
spectrum  exhibited  identifiable  SOj  spectral  absorption  features.  This  demonstrated  that 
the  sensitivity  of  NUVIS  was  insufficient  to  detect  SOj  at  column  abundances  on  the 
order  of  approximately  10^°m'^ 
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Relative  Absorption  ir)  Digital  Number 


Experiment  #1 : 
Sample  Spectra 


igure  4.12.  Sample  spectra  of  experiment  #1  including  measurements  from 
one  evacuated  cell  and  two  SO2  test  cells. 


0.18  • 
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0.14  i 
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Relative  Absorption  Spectra 

(Experiment  #1 :  4"  Test  Cell) 
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580 
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Spectral  Sample  Number 

F igure  4.13.  Relative  absorption  spectra  of  experiment  #  1 . 


E.  GENERATING  EXPERIMENTAL  CURVES  OF  GROWTH 


1.  Approach 

A  second  set  of  experiments  was  undertaken  with  higher  concentrations  of 
SO2.  Various  concentrations  were  tried  until  unambiguous  detection  of  SO2  spectral 
features  was  assured.  Given  the  limitations  of  the  pressure  gage,  described  above,  and  the 
availability  of  only  0.1 1%  and  10%  SO2  mixtures,  the  range  of  reliable  gas  abundances  in 
the  cells  was  somewhat  restricted.  But,  when  data  from  these  abundances  were  combined 
with  those  from  earlier  measurements  using  the  1.55  m  cell,  enough  data  points  for  the 
generation  of  calibration  curves  of  growth  were  achieved.  The  details  are  summarized  in 
Table  4.4.  The  second  experiment  was  initially  designed  to  increase  test  cell 
concentration  in  an  attempt  to  see  the  first  signs  of  absorption  line  formation.  Increasing 
cell  concentration  meant  utilizing  the  10%  SO2  as  the  test  source.  Unfortunately,  the 
restrictions  imposed  by  the  Matheson  gage  setup  (Figure  4.4)  limited  the  four  inch  test 
cell  to  be  backfilled  with  10%  SO2  only  as  low  as  50  mm  Hg.  This  was  the  lowest 
achievable  concentration  with  the  existing  pressure  regulator,  gage  and  valve  manifold. 
Details  of  all  cell  parameters  and  concentrations  utilized  in  generating  the  curves  of 
growth  may  be  found  in  Table  4.4. 


Sample 

Number 

Cell 

Source 

(“/oSOj) 

Pressure  (mm  Hg) 

Path  length  (m) 

lA 

4” 

50 

.0325 

2B 

1.55m 

0.11 

100 

1.55 

3A 

4” 

10 

70 

.0325 

4A 

4” 

10 

100 

.0325 

5B 

1.55m 

0.11 

400 

1.55 

6B 

1.55m 

0.11 

760 

1.55 

Table  4.4.  Cell  parameters  utilized  in  generating  the  curves  of  growth. 
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2.  Experimental  Results 

The  resultant  raw  spectra  of  experiment  #2  are  illustrated  in  Figure  4.14.  When  these 
SO2  spectra  are  ratioed  with  their  corresponding  evacuated  cell  spectra,  the  ratio 
spectra  of  Figure  4.15  are  generated.  These  spectra  are  plotted  along  with  three 
additional  absorption  spectra  (calculated  from  the  Cleary  data  sets)  in  Figure  4.16.  The 
data  sets  are  marked  A  (new)  and  B  (old)  data  respectively.  They  have  been  adjusted 


Raw  Spectra 

(Experiment  #2) 


Spectral  Sample  Number 

Figure  4.14.  Raw  spectra  from  experiment  #2,  including  measurements 
from  one  evacuated  and  three  SO2  test  cells. 


56 


Relative  Absorption 


Relative  Absorption  Spectrum 

(Ratio  Spectra  from  Experiment  #2) 


500  520  540  560  580  600  620  640 


Spectral  Sample  Number 

Figure  4. 15,  Relative  absorption  spectra  from  experiment  #2. 

Relative  Absorption  Spectra 
(Ratio  Spectra  from  Combined  Data  Sets) 


350  400  450  500  550  600 

Spectral  Sample  Number 

Figure  4. 16.  Relative  absorption  spectra  generated  from  combined  data 
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to  the  same  baseline  to  compensate  for  instrumental  gain  drift  effects  and  lamp 
fluctuations.  It  was  also  discovered  that  the  new  data  set’s  spectral  features  were  shifted 
by  five  spectral  sample  numbers  from  the  earlier  data.  This  was  not  too  surprising,  since 
the  instrument  had  now  completed  four  major  field  deployments  and  had 
experienced  great  deal  of  handling  during  shipment.  To  compensate,  bands  in  the  new 
spectra  (Alpha  data  sets)  were  unifonnly  adjusted  by  five  spexels  to  match  the 
corresponding  absorption  lines  of  the  old  spectra  (Beta  data  sets).  Figure  4.16  displays 
three  of  the  new  data  sets  prior  to  adjustment  (one  is  annotated  with  an  arrow)  to 
demonstrate  this  effect.  The  adjusted  spectra  are  illustrated  in  Figure  4. 17. 

Adjusted  Relative  Absorption  Spectra 


Figure  4. 17.  Adjusted  relative  absorption  spectra.  New  data  sets  adjusted  to 
account  for  wavelength  shifts  in  the  instrument.  Please  refer  back  to  Table  4.4  for 
details  on  cell  concentrations  1 A  through  6B. 
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3.  Analysis  of  Experimental  Data 

The  IDL  proceedure  GAUSSFIT  was  utilized  to  fit  individual  Gaussian  profiles  to 
each  of  four  absortion  lines  (centered  about  spectral  sample  numbers  535,  554,  572  and 
590)  in  each  of  the  six  data  sets.  The  resultant  fits  are  illustrated  in  Figures  4.18  through 
4.23.  Their  equivalent  widths  were  calculated  by  integrating  each  Gaussian  (24  in  all)  in 
accordance  with  Equation  2.9.  Table  4.5  provides  details  for  each  of  these  24  absorption 
lines.  For  each  absorption  line,  six  different  equivalent  width  measurements  were 
recorded,  one  for  each  SO2  concentration.  Equivalent  widths  were  then  plotted  against 
cell  concentration  and  resulted  in  the  generation  of  four  curves  of  growth.  A  least  squares 
nonlinear  fit  has  been  performed  to  each  of  the  four  sets  of  measurements  (at  spectral 
sample  numbers  535,  554, 472,  and  590,  respectively)  using  Mathematica  to  produce  four 
modeled  calibration  curves.  A  summary  of  the  resultant  equations  for  each  of  the  four 
models  is  provided  in  Table  4.6  and  curves  of  gro-wth  are  plotted  in  log-log  scale  format 


in  Figures  4.24. 


SO2  Column  Abundance 
(#molecules/m^)  xlO^l 

5.33 

5.61 

7.48 

10.7 

22.5 

42.4 

Nominal 

A,: 

Equivalent  Width  of 
Spectral  Feature  in  nm  at: 

0.382 

0.561 

0.761 

1.12 

1.64 

2.57 

306.5 

nm 

Sample  #590 

308.8 

nm 

Sample  #572 

0.375 

0.396 

0.576 

0.673 

1.06 

1.40 

310.9 

nm 

Sample  #554 

0.250 

0.303 

0.422 

0.739 

0.902 

1.68 

313 

nm 

Sample  #535 

0.211 

0.281 

0.384 

0.505 

1.05 

1.39 

Table  4.5.  Line  equivalent  widths  for  measued  absorption 

column  abundances. 


lines  as  a  fimction  of  SO2 
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Rol»il«o  Absorption  RelaUve  Absorption  Relative  Absorption 


Gaussian  Fits 

(Spectral  Sample  Numbers  =  Wavelengths) 


Spectral  Sample  Number  Spectral  Sample  Number 


Figure  4. 18.  Concentration  1 A  Figure  4. 19.  Concentration  2B 


Figure  4.20.  Concentration  3A  Figure  4.2 1 .  Concentration  4A 


0 

eCD5Z)5O5e0SD600  6Q06O 


Figure  4.22.  Concentration  5B  Figure  4.23.  Concentration  6B 
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Sample  Number 
of  Spectral 
Feature 

Equation  for  Modeled  Curves  of  Growth 

535 

1F= -28.234 +  0.5685  In(^) 

554 

-31.122 +  0.6268  In(^) 

572 

-23.975  +  0.4867  In  (^) 

590 

1F= -48.372 +  0.9754  In  ($•) 

Ta 

Die  4.6.  Equations  for  modeled  curves  of  growth. 

(I^is  the  equivalent  width  in  nm  and  ^is  the  column  abundance  in  m'^) 


5x10*'  1x10**  2x10**  5x10** 

SOj  Column  Abundance  (m"^) 

Figure  4.24.  Laboratory  generated  curves  of  grovvth  for  SO2 
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V.  APPLICATION  OF  LABORATORY  CALIBRATIONS  TO  FIELD  DATA 

As  stated  earlier,  the  motivation  for  this  thesis  research  has  been  the  need  to 
characterize  and  to  quantify  NUVIS’s  capabilities  as  a  remote  sensing  spectrochemical 
analysis  instrument.  The  completion  of  the  laboratory  calibration  procedures  described  in 
Chapter  IV  now  permit  us  to  apply  these  quatitative  calibrations  to  real  field  data. 

A.  FIELD  EXPERIMENTS 

Data  described  in  this  chapter  resulted  from  three  separate  field  experiments  with 
NUVIS.  Two  field  collections  were  conducted  in  1998,  just  prior  to  the  inception  of  this 
thesis  project.  The  third  field  collection  was  conducted  in  April  of  1999  in  direct  support 
of  this  thesis  research.  Each  of  these  experiments  was  conducted  on  site  at  a  coal-firing 
power  plant.  The  plant  provided  three  operational  smokestacks  that  were  available  for 
remote  measurement  through  clear  desert  air. 

During  the  1998  collections,  only  two  of  the  three  operational  stacks  had 
operational  SO2  scrubbers.  Thus,  the  third  stack  provided  a  much  higher  SO2  plume 
signature  than  the  other  stacks.  Direct  in  situ  measurements  provided  by  plant  managers 
indicated  that  the  clean  stacks  emitted  SO2  at  levels  in  the  tens  of  ppm  (approximately  40 
ppm).  Similar  measurements  indicated  that  the  unscrubbed  stack  operated  with  SO2 
emissions  in  the  hundreds  of  ppm  (approximately  100  to  400  ppm).  Photos  of  the  power 
plant  may  be  found  in  Figure  5.1.  There  are  actually  six  stacks  in  the  image  scene.  Three 
of  the  six  were  recently  built  because  of  mandated  reductions  in  pollutant  emissions. 
(Cleary,  1998) 


63 


Figure  5.1.  Coal  burning  power  plant.  Photos  were  taken  during  1998  (left)  and 
1999  (right)  collections.  Only  two  of  the  three  scrubbers  were  in  operation  in  1998 
but  in  1999  scrubbers  were  on-line  in  all  three  stacks. 

1998  photo  courtesy  of  Cleary,  1998. 


During  the  most  recent  on  site  field  measurements,  in  the  spring  of  1999,  the  power  plant 
had  just  brought  the  final  set  of  scrubbers  on-line.  This  meant  that  all  three  stacks 
emitted  SO2  at  or  below  40  ppm  concentration. 

The  acquisition  of  experimental  field  data  followed  a  procedure  much  like  that 
used  in  the  laboratory,  but  with  the  following  modifications:  1)  NUVIS  was  operated  in 
its  full  hyperspectral  imaging  mode,  providing  spectral  information  for  each  2D  image 
pixel.  2)  The  source  of  backgroimd  continuum  radiation,  against  which  SO2  absorption 
was  measured,  was  natural  sky  scattered  UV  solar  radiation.  Hence,  each  spectrum  had 
to  be  ratioed  with  a  clear  sky  spectrum  in  order  to  quantify  SO2  absorption.  Figure  5.2 
shows  a  representative  sky  spectrum  recorder  in  the  field. 

Figure  5.3  shows  a  corresponding  spectrum  for  a  region  of  an  image  including  the 
effluent  plume  from  a  smokestack.  It  is  evident  that  absorption  features  exist  that  are  not 
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present  in  the  clear  sky  spectrum.  This  is  particularly  obvious  in  Figure  5.4,  where  both 
spectra  are  plotted  jointly.  Figure  5.5  shows  the  ratio  of  the  plume  and  clear  sky  spectra. 
The  typical  SOj  spectral  signatures  are  evident  at  the  high  spectral  sample  side  of  the 
plot.  By  fitting  profiles  to  these  signatures,  calculating  each  profile’s  equivalent  width, 
and  comparing  those  with  the  calibration  curves  of  grov^h,  effective  column  abtmdances 
of  SO2  can,  in  principle,  be  determined. 


Sky  Data 


Spectral  Sample  Number  (Spexel) 


Figure  5.2.  Sky  spectrum  recorded  by  NUVIS. 
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Plume  Data 


Ratio  of  Plume  to  Sky  Spectrum 


Figure  5.5.  Plume-to-sky  ratio  spectrum.  Values  range  from  0-1.0 
corresponding  to  0-100%  absorption. 


B.  RESULTS 

In  order  to  test  the  calibration  procedure,  representative  sky  and  plume  regions  of 
interest  (ROIs)  were  selected  for  the  “dirty”  stack  from  one  of  the  1998  data  sets.  This 
particular  set  was  obtained  at  a  range  of  approximately  0.2  km  from  the  stack.  Figure  5.6 
illustrates  the  regions  of  interest  selected  for  analysis.  Figure  5.7  illustrates  the  mean 
spectra  for  each  region,  and  Figure  5.8  displays  the  resultant  ratio  absorption  spectrum  for 
the  plume,  which  has  been  labeled  as  stack  #1.  These  results  validated  the  claim  that 
NUVIS  had  the  capability  of  detecting  sulfur  dioxide  plumes  at  short  range.  The  image 
of  Figure  5.7  is  produced  through  the  Regions  of  Interest  algorithms  of  the  Environment 
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for  Visualizing  Images  (ENVI)  software.  Additional  ENVI  analysis  is  further  detailed 
in  Appendix  C. 


Stack  #1  with  ROI 
(at  0.2km) 

■  Sky  (1000  pts) 
g  Plume  (1000  pts) 


Figure  5.6.  Regions  of 
Interest  (ROI’s)  used  to 
extract  sky  and  plume 


spectra. 


ROI  Mean  Spectra 


Figure  5.7.  ROI  mean  spectra 
(see  ROI’s  of  Figure  5.6) 


Relative  Plume  Absorption 


Figure  5.8.  Relative  plume  absorption, 
(see  Figure  5.6) 
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Applying  the  same  Gaussian  curve  fitting  techniques  to  this  spectrum  as  were  used  for 
the  lab  data  resulted  in  the  curve  fit  illustrated  in  Figure  5.9.  Equivalent  width 
measurements  were  made  for  each  of  the  absorption  lines  centered  at  spectral  sample 
numbers  535,  554  and  572,  respectively. 


Stack  #1  Gaussian  Fit  (at  535, 554, 572) 


This  resulted  in  equivalent  widths  of  1.80,  1.40,  and  2.03  nm,  respectively.  The 
laboratory  curves  of  growth  were  then  used  to  deduce  the  corresponding  plume  SO2 
column  abundances.  Results  from  this  application  of  the  calibration  curves  to  the  field 
data  are  summarized  in  Table  5.1. 
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Spectral  Sample 
Number 

Equivalent  Width;  W 
(nm) 

Column  Abundance:  c, 
(m-^) 

535 

1.80 

8.88  X  10^^ 

554 

1.40 

3.42  X  10^^ 

572 

2.03 

Table  5.1.  Resultant  column  abundances  of  stack  #  1  j 

fom  three  spectral  features. 

These  are  the  effective  column  abundances  for  the  plume,  assuming  that  it  is  a 
layer  of  gas  of  uniform  number  density  n  and  thickness  L,  in  which  case  the  column 
abundance  would  be  simply 

g=nL.  (5-1) 

However,  the  actual  plume  geometry  is  approximately  cylindrical,  rather  than  layer-like. 
So  the  actual  effective  absorption  path  length,  assuming  a  uniform  effluent  density  across 
the  plume  of  diameter  D,  would  be  the  mean  thickness  of  a  cylinder  of  that  diameter, 
as  illustrated  in  Figure  5.10  and  in  Equation  5-2  where  the  effective  path  length  is  L^, 


STACK  PATH  LENGTH 


NUVIS 


> 


Figure  5.10.  Plume  geometry. 
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(5-2) 


-0/2 


Therefore,  for  a  measured  column  abundance  g  across  a  plume  of  diameter  D,  the 
number  density  of  SO2  in  the  plume  will  be 


(5-3) 


Of  course,  this  assumes  that  the  spectrum  was  measured  along  a  direction  perpendicular 
to  the  cylinder  axis  of  the  vertical  plume  (in  other  words,  in  a  horizontal  direction).  If  the 
measurement  is  perfomied  at  close  range,  sighting  at  an  elevation  angle  0  with  respect  to 
the  horizontal,  then  be  projected  along  a  slant  path,  giving 


kD 

4cos^’ 


(5-4) 


from  which  it  follows  that  the  measured  SO2  number  density  will  be 


4$'C0S^  (5-5) 

kD 
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Spectral  Sample  Number 

SO2  Number  Density  (m'^) 

535 

6.75  X  10^* 

554 

572 

1.22x10^^ 

Table  5.2.  SO2  number  densities  with  geometric  scaling. 


The  spectrochemical  analysis  just  discussed  yielded  the  plume  SO2  abundance  in 
the  form  of.  a  number  density  (i.e.  number  of  molecules  per  cubic  meter  of  plume  gas). 
However,  in  situ  pollution  monitoring  devices  in  the  stacks  express  SO2  concentrations  in 
parts  per  million  (ppm)  of  stack  effluent,  a  standard  unit  of  measure  in  the  fields  of 
chemical  and  environmental  engineering,  called  the  gas  mixing  ratio.  So,  we  must 
convert  our  spectrochemical  measurements  to  those  units  if  a  meaningful  comparison  is 
to  be  made. 

Images  indicate  that  cylindrical  plumes  rise  a  considerable  distance  above  the 
tops  of  the  stacks  without  changing  their  overall  shapes  or  diameters.  This  implies  that 
the  gas  pressure  in  the  flow  must  be  approximately  isobaric  with  the  surrounding  cool  air. 
This  means  that,  for  the  comparitively  short  section  of  plume  observed  immediately 
above  the  flues,  Bernoulli’s  equation  (Halliday  and  Resnick,  1997)  demands  that 


•  2 

Pair~  P plume  0-5  PpUime  ^  plume : 


(5-6) 


where  Pair  is  the  hydrostatic  pressure  of  the  surrounding  cool  air,  PpUime  is  the 
corresponding  hydrostatic  pressure  within  the  warm  plume,  PpUime  is  the  mass  density  of 
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the  plume  gas,  and  Vpiume  is  the  plume  gas  flow  velocity  with  respect  to  the  surrounding 
air. 

The  actual  situation  is,  to  a  good  approximation,  somewhat  simpler  than 
Bernoulli’s  equation  would  imply.  Since 

Pair  =  l.OxlO'Pa,  pplume=Pair  =  l-3kg  m'^  and  VpUrrne  =  10ms'\  (5-7) 

the  magnitude  of  the  hydrodynamic  term  in  Bemouilli’s  equation  is  roughly 

0-5  Pplume^  plum^  65Pa.  (3"8) 

This  is  much  smaller  that  Pair  and  so  may  be  neglected  in  the  analysis.  So,  we  have  a 
quasi-hydrostatic  system  in  which  Ppiume  equals  Pair^o  a  good  approximation.  Ppiu^e  must 
be  the  sum  of  the  partial  pressures  Pso2  and  Pother  of  the  SO2  and  other  effluent 
constituents,  respectively.  Assuming  that  the  plume  is  well  mixed  and  is  in 
thermodynamic  equilibrium  at  absolute  temperature  Tphme,  the  ideal  gas  law  (Halliday 
and  Resnick,  1997),  when  combined  with  the  above  results,  requires  that 

nso2k  1'plume'^  Mother  k  Ppiume  ^air  k  Pair  5  (5-9) 

where  nso2,  ngther,  Hair  are  the  number  densities  of  SO2,  other  plume  constituents,  and  the 
surrounding  air,  respectively,  k  is  Boltzman’s  constant  (1.38  x  10'^^  J/K),  and  Pair  is  the 
absolute  temperature  of  the  surrounding  air.  Therefore,  it  follows  that 

Mother  (f^air  Pair  I  Ppiume)  -  ns02>  (5-10) 

and  for  each  spectrochemical  measurement  of  nso2  ,  this  formula  may  be  used  to 
calculate  Hother  ■  With  that  number  in  hand,  the  measured  plume  gas  mixing  ratio 

nso2/nother  (5-11) 
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may  be  estimated  and  compared  with  the  in  situ  SO2  measurements.  Here  are  the  results: 


Spectral  Sample 
Number 

nS02 /Mother 

Mixing  Ratio  (ppm) 

535 

3.43  xio"* 

343 

554 

1.32x10“^ 

132 

572 

6.20  xlO-^ 

620 

Assumed  Values; 

tlair 

2.69  xlO^"  m‘^ 

'^plume 

410K 

ttair’I'aii^ plume 

1.968x10^^  m'^ 

Table  5.3.  Mixing  ratio  results  from  observations  of  stack#l  at  0.2  km. 


So  the  average  SO2  mixing  ratio  is  365±  200ppm,  where  the  uncertainty  is  simply 
the  standard  deviation  about  the  average.  This  agrees  very  well  with  the  in  situ  direct 
measurements  which  are  near  400  ppm.  Because  we  have  only  three  measurements  to 
work  with,  the  statistics  are  not  overwhelmingly  compelling,  but  the  consistency  between 
the  two  independent  measurements  is  still  quite  promising. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

This  thesis  research  demonstrated  conclusively  that  NUVIS  is  capable  of 

detecting  and  quantifying  SO2  emissions  in  the  form  of  effluent  smokestack  plumes.  This 
effort  entailed  two  separate  phases  of  research.  First,  laboratory  comparison  UV  spectra 
of  SO2  were  recorded  at  different  concentrations.  These  spectra  were  used  to  calculate 
least-squares  fitted  curves  of  growth  for  four  different  SO2  spectral  features.  The  curves 
of  growth  can  now  be  used  for  any  subsequent  spectrochemical  measurements  of  SO2 
column  abundances  by  NUVIS.  Second,  the  results  of  the  laboratory  studies  were 
employed  to  analyze  real  field  data  from  a  coal-burning  power  plant.  Direct  in  situ 
chemical  measurements  of  the  SO2  in  one  of  the  plant’s  effluent  plumes  was  known  to  be 
400  ppm.  Detailed  analysis  of  NUVIS  data  from  the  plume  yielded  a 'mean  plume  SOj 
mixing  ratio  of  365  ±  200  ppm,  in  excellent  agreement  with  the  in  situ  value. 

Based  upon  an  assessment  of  the  signal-to-noise  ratios  in  representative  field  data, 
and  the  laboratory  observation  that  SO2  column  abundances  must  be  >  10^°  m'^  for 
NUVIS  to  detect  the  SO2  signatures  that  we  have  analyzed,  we  can  estimate  the  practical 
lower  mixing  ratio  limit  for  SO2  detection  in  the  field.  We  estimate  that  it  is 
approximately  70  ppm.  This  indicates  that  NUVIS  should  not  have  detected  quantifiable 
SO2  emissions  during  the  1999  observations  of  the  power  plant,  when  all  of  the  plant’s 
plumes  were  chemically  scrubbed  down  to  SO2  mixing  ratios  on  the  order  of  40  ppm. 
That  is  indeed  the  case,  further  demonstrating  the  consistency  of  the  experimental  and 
analytical  procedures  developed  during  this  thesis  research  project. 
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A  number  of  steps  could  be  taken  to  enhance  the  calibration  and  analysis 
techniques  outlined  in  this  thesis  and  to  improve  the  overall  sensitivity  of  the  UV 
hyperspectral  remote  sensing  instrumentation.  1)  A  ne'w  wavelength  calibration  should 
be  performed.  2)  A  new  set  of  calibration  curves  should  be  developed,  involving  a  wider 
range  of  SOj  concentrations.  3)  Repetitive  measurements  should  be  made  at  each 
location  for  each  scene  when  in  the  field  in  order  to  check  for  consistency  and  to  improve 
upon  the  signal-to-noise  in  the  data.  4)  Further  develop  and  extend  the  analysis  processes 
presented  in  Chapters  IV  and  V  to  include  additional  spectral  features,  thereby  providing 
a  larger  statistical  set  of  measurements.  5)  NUVIS  should  be  augmented  with  peripheral 
devices  to  measure  all  applicable  environmental  factors  (i.e.  barometric  pressure, 
temperature,  GPS  generated  position)  in  conjimction  with  each  field  measurement. 

As  a  follow-on  to  this  thesis  work,  more  rigorous  analysis  should  be  conducted  on 
the  existing  field  data  utilizing  the  algorithms  provided  by  the  ENVI  software  to  more 
reliably  quantify  the  column  densities  through  a  more  refined  application  of  this 
calibration.  This  process  should  be  further  developed  through  continued  field  and 
laboratory  experimentation  and  should  also  involve  corrections  for  UV  foreground 
scattering  and  the  contributions  due  to  other  chemical  effluents. 

This  study  has  provided  an  account  of  NUVIS  progression  to  date.  It  also 
demonstrated  that  sulfur  dioxide  does  exhibit  a  characteristic  ultraviolet  absorption 
signature.  It  also  indicated  that  the  existing  instrument  was  unable  to  detect  SOj 
signatures  at  concentrations  of  tens  of  ppm  at  long  ranges  (see  Appendix  C). 
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The  sensor’s  resolving  povfer  is  limited  by  the  eight  bit  processing  capability  of 
the  detector,  a  Pulnix  CCD  camera.  Instrument  anomalies  (also  briefly  detailed  in 
Appendix  C)  can  most  likely  be  attributed  to  optical  components  such  as  the  scanning 
mirror  and  diffraction  grating.  The  anomalies  experienced  with  the  processing  PC  can 
most  likely  be  attributed  to  the  instrument’s  software.  In  either  case,  further  testing  is 
required  to  fully  determine  the  extent  of  such  anomalies. 

Further  improvements  to  the  existing  instrument  could  be  made  to  correct  many 
of  these  problems  but  a  more  beneficial  approach  involves  construction  of  a  new 
instrument,  incorporating  the  lessons  learned  from  NUVIS.  It  should  have  higher 
resolving  power  and  greater  sensitivity  utilizing  a  digital  twelve/sixteen  bit  CCD  camera, 
more  robust  processing  power  and  an  improved  electronic  design  that  minimizes  spurious 
operational  anomalies  and  glitches.  We  recommend  that  further  improvements  to  NUVIS 
be  made,  and  that  a  new,  more  capable  instrument  be  developed  and  tested  with  the 
ultimate  goal  of  someday  fielding  a  similar  airborne  UV  spectral  imager. 
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APPENDIX  A.  NUVIS  FIELD  CHECKLIST 

Field  experiments  involve  a  great  deal  of  logistics  to  accomplish  a  successful  mission. 
Ensuring  that  NUVIS  and  all  its  support  components  arrive  at  the  site  is  only  one  of  many 
concerns  during  such  operations.  The  following  checklist  is  provided  as  a  material  guide 
for  NUVIS  field  operation.  The  list  is  by  no  means  all-inclusive  and,  depending  on  the 
mission  at  hand  may  require  modifications.  Five  shipping  crates  (including  pre¬ 
designated  packing/shipping  material:  i.e.  foam,  bubble  wrap  etc...)  are  dedicated  to  the 
field  deployment  of  NUVIS.  Each  shipping  case  will  first  be  listed  followed  by  the 
recommended  contents  of  each  crate; 


•  CASE  #1: 

1)  CPU  (IDL  HASP  may  remain  in  place). 

2)  Monitor  (monitor  sits  firmly  against  right  wall  of  the  shipping  crate  with  foam  installed). 

3)  Max-4 10&  JS  auxiliary  connector  terminator. 

[Notes:  Foam  has  been  specifically  cut  to  specification  for  packing  purposes.  Foam  pieces  are 
labeled  and  fit  nicely  into  place  to  allow  shipping  of  the  CPU,  monitor  and  Max-410  remain  in 
one  crate.  It  is  recommended  that  only  the  foam  be  removed  and  that  these  components  remain  in 
the  crate  during  field  operation.] 


•  CASE  #2 

1)  NUVIS  tripod  (with  wrench). 


•  CASE  #3 

1)  Spare  PC. 
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•  CASE  #4 


1)  25  foot  NUVIS  umbilical  cable. 

2)  8  foot  NUVIS  umbilical  cable. 

3)  Extension  cords  (quantity  is  mission  dependent). 

4)  Category  5  cables,  blue  in  color  (2  jumpers,  2  network). 

5)  Monitor  cable. 

6)  COMl  to  Max-410  interface  cable. 

7)  UPS. 

8)  Power  strip. 

9)  Surge  suppressor. 

10)  Mouse. 

1 1)  Tool  kit  (screwdrivers,  utility  knife,  wire  cutters,  pliers,  Allen  wrenches,  etc...). 

12)  Miscellaneous  items  (paper  towels,  rope,  lens  paper,  gloves,  plastic,  etc...). 


•  CASE  #5 

1)  Spare  stepper  motor. 

2)  Spare  Pulnix  camera  &  Lens. 

3)  Spare  AZ  encoder. 

4)  Spare  computer  keyboard. 

5)  Spare  mouse. 

6)  Spare  JAZ  drive. 

7)  Spare  IDL  HASP. 

8)  JAZ  Cartridges. 

9)  Diagnostic  tools.  Fluke,  meter  and  oscilliscope. 

10)  Repair  items  (tape,  cable  ties,  soldering  iron,  solder,  sponge,  spare  wire,  etc...). 

1 1)  Documentation  (software  documentation,  Keithly,  US  Digital  and  Max-410  docs.). 

12)  Backup  discs  (operating  system,  software,  hard  drive,  etc...). 

13)  Gas  test  cells. 

,  14)  Spare  camera. 


80 


APPENDIX  B.  NUVIS  OPERATIONAL  CHECKLIST 


SETTING  UP  THE  SYSTEM 

#1  Remove  tripod  from  case  and  set  in  place.  Use  key  (stored  at  the  base  of  the 
tripod)  to  move  instrument  locking  mechanism  from  its  slot  by  unlocking  the  setscrew. 

#2  Remove  instrument  from  case  and  slide  instrument  into  place  on  the  tripod 
moimting  plate.  Once  instrument  is  in  place,  reset  locking  mechanism  and  retighten  set¬ 
screw.  (Key  may  be  restowed.) 

#3  Remove  protective  cover  from  optical  sight.  It  is  recommended  that  the  window 
cover  not  be  removed  until  system  is  fully  operational.  This  will  provide  maximum 
protection  to  the  filter  window. 

#4  Ensure  computer  system  components  are  accessible.  Currently  it  is  recommended 
that  the  computer  components  and  stepper  motor  remain  in  the  carry  case  with  only  the 
two  access  covers  removed.  This  provides  maximum  protection  and  ease  of  mobility 
while  all  components  remain  readily  accessible. 

#5  Connect  the  UPS  to  the  appropriate  source  of  power  (115  VAC  outlet  or  field 
generator).  Connect  the  computer,  monitor,  Ethernet  hub  and  stepper  motor  controller 
power  cables  to  each  component  and  to  the  UPS.  Ensure  that  the  voltage  switch  on  the 
computer  is  set  to  1 15  volts. 

#6  Connect  monitor  cable,  mouse'and  keyboard  to  the  PC. 

#7  Now  connect  the  remaining  cables  including  NUVIS  umbilical  cables  and  blue 
category  five  cables.  Connections  are  to  be  made  to  the  varios  cable  connectors  on  the 
NUVIS  instrument  as  well  as  to  the  AMS  410  and  rear  panel  of  the  PC.  Figure  B.l 
illustrates  these  connections. 


Figure  B.l.  NUVIS  umbilical,  instrument,  PC  and  AMS-41U  cable  connections. 
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GETTING  UP  AND  RUNNING 


(You  should  now  be  ready  to  power  up  the  system.) 

#1  Turn  on  the  monitor  and  PC.  (To  turn  on  PC  ensure  switch  on  back  panel  is  set  to 
the  on  position  and  then  push,  turn  and  release  blue  switch  on  the  front  panel  of  the  PC.) 

#2  Wait  for  the  BIOS  display  and  initial  indication  that  the  operating  system 
(Windows  95)  is  functioning.  This  will  take  some  time. 

#3  Once  this  is  complete,  the  standard  ‘Enter  Network  Password’  screen  will  appear. 
If  networking  is  desired,  enter  appropriate  password,  otherwise  press  ‘escape.’  In  either 
case,  the  standard  Windows  screen  with  short-cut  icons  will  next  appear. 

DO  NOT  LAUNCH  THE  NUVIS.EXE  PROGRAM 
UNTIL  COMPLETING  STEPS  SIX  THROUGH  EIGHT: 

#6  Turn  on  the  stepper  motor  controller. 

#7  Once  this  is  accomplished  and  you  are  ready  to  take  data,  remove  the  protective 
cover  from  the  NUVIS  window.  (If  the  4inch  test  cell  is  required,  install  the  mounting 
studs,  install  the  cell,  and  secure  the  cell  in  place  with  the  stud  nuts.)  Align  the 
instrument  as  required. 

#8  The  NUVIS.EXE  program  is  memory  intensive.  Close  all  other  applications 
before  launching  this  program.  It  is  actually  recommended  that  NO  OTHER 
APPLICATIONS  BE  OPENED  PRIOR  TO  LAUNCHING  NUVIS.EXE.  If  a  number  of 
applications  had  been  previously  opened,  it  might  be  wise  to  shut  down  the  PC  and 
restart.  This  will  dedicate  the  PC  to  NUVIS  operation  only.  Once  this  is  complete,  you 
are  now  ready  for  NUVIS  operation. 


INSTRUMENT  OPERATION  WITH  NUVIS.EXE 
#1  Launch  the  NUVIS.EXE  program. 

#2  Upon  execution  of  the  program,  four  windows  will  appear  on  the  monitor  prior  to 
the  presentation  of  the  actual  operational  window.  These  windows  are  detailed  in  Table 
B.l. 
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SEQUENCE 

WINDOW  STATEMENT 

DESCRIPTION/ 

REQUIRED  ACTION 

1) 

Is  Max-410  Ready? 

This  is  a  reminder.  As  long  as  Max-410  is  on,  enter 

OK. 

2) 

Mirror  is  aligned. 

Once  previous  step  was  completed,  the  Max-410 
ensured  communications  with  the  instrument  and 
centered  the  mirror  to  position  zero.  This  statement 
confirms  that  action.  Enter  OK. 

3) 

Blank 

Blank  Window.  Enter  OK 

4) 

Based  on  memory,  #  frames 
is  XXX.  (i.e.  613) 

This  indicates  the  number  of  images  that  may  be  stored 
in  RAM.  If  this  number  is  below  400,  another 
application  was  most  likely  previously  opened.  Before 
proceeding,  it  is  best  to  restart  the  sequence  to  ensure 
enough  memory  is  available.  Otherwise,  Enter  OK. 

Table  E 

>.  1 .  Window  sequence  upon  execution  of  the  NUVIS.EXE  program. 

#3  Once  the  above  sequence  is  completed,  the  operational  window  of  Figure  B.2  is 
obtained  and  preparation  for  an  actual  data  collection  is  nearly  complete. 


Figure  B.2.  NUVIS  operational  Window. 
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(If  any  problems  are  encountered  prior  to  this  point,  it  is  best  to  exit  the  program  and 
restart  the  process.  If  memory  was  initially  low,  this  may  require  a  complete  system 
shutdown  and  restart.  If  the  program  screen  locks  up  during  portion  of  the  startup 
sequence,  the  'ctrl-alt-del'  command  to  end  the  task  may  not  be  available.  If  this  occurs, 
do  not  attempt  to  use  this  command  a  second  time  to  restart  your  computer.  This  method 
of  rebooting  (once  the  NUVIS.EXE  program  is  executed)  may  cause  damage  to  certain 
Windows  operating  files.  As  a  last  resort,  if  the  program  locks  up,  try  to  turn  the  system 
power  off  manually  and  restart  the  system.  Normal  system  shutdown  is  detailed  in  the 
final  steps  of  this  procedure.  Likewise,  when  exiting  the  NUVIS.EXE  program,  it  is  best 
to  perform  a  proper  system  shutdown  every  time.  This  requires  using  the  EXIT  button 
appearing  in  the  operational  control  window. 

#4  Set  the  Photomultiplier  Gain.  It  is  best  to  start  at  a  lower  value  (5.0  or  below)  to 
begin.  Gains  may  be  set  from  0  to  10  by  keyboard  entry  or  by  using  the  arrows  and 
mouse. 

#5  Assuming  that  the  instrument  is  pointed  at  the  appropriate  target,  a  quick-look  at 
the  sensor’s  active  digital  signal  may  be  obtained  by  activating  the  VIDEO  button.  The 
brighter  pixels  correspond  to  higher  digital  signals.  Keeping  in  mind  that  the  detector 
should  not  be  saturated,  the  image  intensifier  should  be  adjusted  until  an  appropriate 
signal  is  displayed.  The  image  intensifier  must  also  first  be  turned  ON  prior  to  selecting 
the  VIDEO  button  or  either  of  the  two  modes  described  in  the  steps  that  follow.  Also, 
keep  in  mind  that  the  screen  display  does  not  represent  the  image  scene  directly  as 
viewed  from  behind  the  instrument.  The  screen  display  is  actually  inverted  from  the 
scene  itself 


#6  Prior  to  collecting  data,  ensure  that  there  is  enough  memory  available  on  the 
hardrive  and/or  ensure  that  a  JAZ  cartridge  is  installed  for  data  collection. 


#7  Mirror  Position  (center  position)  may  be  set  through  the  active  window.  This 
number  should  initially  be  zero  since  the  mirror  was  centered  during  the  startup  sequence. 

#8  Two  operational  modes  exist  with  NUVIS;  (1)  SNAP  and  (2)  SCAN.  The 
following  information  should  be  reviewed  before  selecting  either  one  of  these  modes: 

(1)  SNAP: 

Images/Snap.  This  setting  determines  the  successive  number  of  SNAPS  to  be  performed. 
During  a  NUVIS  SNAP,  the  mirror  remains  centered  on  the  position  of  step  #7. 

(2)  SCAN: 
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The  ^tart  and  Stop  settings  only  apply  to  the  SCAN  function.  Setting  mirror  Start  and 
Stop  positions  automatically  resets  the  number  of  steps  through  which  the  mirror  will 
scan.  The  number  of  steps  can  only  be  controlled  by  the  start  and  stop  settings.  A 
maximum  of  351  steps  may  be  utilized  but  steps  in  excess  of  200  may  cause  system 
failures  due  to  memory  allocation  faults. 

The  number  of  Images/Step  may  also  be  selected.  This  defines  the  number  of  SNAPS  to 
be  performed  at  each  repositioning  of  the  mirror 


A  NUVIS  "SNAP"  is  simply  a  640  x  480  byte  array  of  data  achieved  through  the 
instrument  control  program.  Each  individual  SNAP  may  be  written  to  either  the  PC’s 
hard  drive  or  internal  JAZ  drive  and  will  register  as  a  300KB  file  if  only  one  frame  is 
obtained.  One  frame  is  simply  a  two-dimensional  set  of  data  representing  the  spectral 
(wavelength)  and  spatial  (vertical)  components  of  the  representative  scene.  Although  this 
frame  of  data  is  not  really  an  image,  the  operational  NUVIS  window  refers  to  it  as  such. 

Selecting  multiple  images  per  SNAP  (i.e.  100),  will  increase  the  size  of  the  data  file  as 
100  different  640x480  byte  arrays  are  written,  one  by  one,  to  the  either  of  the  available 
drives  as  a  complete  data  file.  This  is  the  preferred  mode  of  operation  during  an 
instrument  calibration  when  the  entire  field  of  view  of  the  instrument  is  uniformly 
focused  on  a  standard  gas  test  cell  and  not  a  true  image  scene.  These  multiple  SNAPS 
may  be  averaged  during  the  analysis  process  to  improve  signal  to  noise.  A  file  containing 
100  “images  per  SNAP”  is  normally  30MB  in  size.  With  these  settings,  a  data  file 
smaller  than  30MB  indicates  a  problem  with  that  particular  collection. 

A  NUVIS  "SCAN"  is  a  bit  different.  It  is  the  preferred  method  when  an  image  scene  is 
being  targeted.  By  selecting  the  mirror  stop  and  start  positions  in  the  instrument  control 
window  (see  Appendix  B)  the  operator  sets  the  number  of  steps  (i.e.  scan  line  numbers) 
through  which  the  mirror  will  scan.  The  mirror,  which  is  first  centered  at  position  zero 
during  startup,  is  moved  to  the  initial  offset  position  (i.e.  start  position)  and  begins  to 
scan.  It  rotates  through  the  specified  number  of  steps  until  reaching  the  magnitude  of  the 
same  offset  position  (i.e.  stop  position)  at  the  end  of  the  scan.  The  number  of  images  per 
step  may  also  be  set  by  the  operator.  Typical  settings  for  a  routine  SCAN  include 
start/stop  values  of  -50  to  +50,  which  corresponds  to  101  different  steps  or  mirror 
positions.  Approximately  300KB  of  data  storage  is  required  for  each  step.  As 
demonstrated  by  each  of  these  operational  modes,  NUVIS  data  files  are  quite  large 
requiring  that  the  CPU  be  fully  dedicated  to  the  reading/writing  of  this  data  once  a  data 
collection  has  been  initiated.  This  is  a  major  restriction  imposed  by  the  current 
NUVIS.EXE  program  that  is  written  in  Visual  Basic. 
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#9  Select  the  appropriate  mode:  SNAP  or  SCAN. 

#10  Once  the  mode  is  selected  and  executed,  the  reading  and  writing  of  sensor  data 
will  take  a  considerable  amount  of  time  (30  to  120  seconds  in  most  cases).  The  program 
will  ask  for  a  file  name  in  which  to  store  the  data.  The  customary  file  extension  used  in 
all  NUVIS  files  is  a  ‘.nuv’ extension.  The  file  will  be  saved  in  JPEG  format. 

#11  It  is  highly  recommended  that  one  operation  be  fiilly  completed  before  the  next 
command  is  issued.  It  is  best  to  be  patient,  as  the  program  has  its  own  inherent 
limitations  and  is  not  as  streamlined  as  many  of  the  commercial  software  applications  on 
the  market.  System  operators  should  ensure  that  only  one  function  at  a  time  is  requested 
to  reduce  the  possibility  of  program  failures. 

#12  As  previously  mentioned,  to  exit  the  system  use  the  EXIT  button. 

The  TEST  ARC  and  LOAD  IMAGE  buttons  are  not  functional. 
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APPENDIX  C:  RANGE  EFFECTS  AND  NUVIS  ANOMALIES 

Data  presented  in  this  section  were  taken  from  the  1998  field  trip  to  the  coal 
burning  power  plant.  Figure  C.l  illustrates  the  scene  in  a  principle  components  format. 
A  principle  components  transform  on  bands  500  to  640  reduces  the  data  to  this  image, 
which  highlights  the  differences  between  the  clean  and  dirty  plumes.  A  brief  analysis  of 
the  plumes  is  done  with  eight  regions  of  interest  at  a  range  of  1 .5  km.  These  regions  are 
illustrated  in  Figure  C.2.  The  ratioed  relative  absorption  spectra  are  displayed  in  Figure 
C.3,  demonstrating  not  only  the  effects  of  range  on  the  detected  spectra  but  also  providing 
a  comparison  of  the  signatures  of  the  scrubbed  and  unscrubbed  stacks  at  this  range. 

A  very  preliminary  analysis  of  these  data  was  conducted  using  a  matched  filter 
technique.  Because  the  absorption  lines  for  the  clean  stacks  were  not  being 
unambiguously  detected  by  the  instrument  at  range,  this  approach  was  tried  in  order  to 
improve  the  detection  of  SO2  signatures.  The  matched  filter  technique  utilized  in  scene 
spectra  for  both  the  sulfur  dioxide  and  the  sky  backgroimd.  The  results  of  this  technique 
are  displayed  in  Figures  C.4  and  C.5.  In  Figure  C.4,  apparent  regions  of  SO2  are 
indicated  by  the  dark  regions  in  the  display.  This  particular  scene  includes  the 
unscrubbed  stack  at  a  range  of  1.5km.  A  similar  matching  technique  resulted  in  pixel 
matching  in  the  modeled  sky  spectrum  as  illustrated  by  the  lighter  regions  of  Figure  C.5. 
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3  Stacks  - 1.5  km 
Principal  Components 

Cyan:  PC1  Red:  PC2 


Figure  C.  1.  Principle  components  image 
taken  from  bands  500-640  of  the  original 
data  set.  Principal  component  1  is  encoded 
as  cyan  (green/blue  guns),  principal 
compent  2  is  encoded  as  red. 


Regions  of  Interest 
at  1 .5  km 


I  Stack  1A 
■  Stack  IB 


M  Stack  1D 

■  Stack  2A 

■  Stack  2B 

■  Stack  3A 

■  Stack  3B 

■  Sky 


Figure  C.2.  Principle  component  1  image 
with  regions  of  interest  superimposed. 

Color  codes  correspond  to  the  mean  spectra 
in  Figure  C.3. 
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Relative  Absorption 


Figure  C.4.  SO2  matched  filter.  Scale  Figure  C.5.  Sky  matched 
indicates  relative  match.  filter. 
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